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^ . We show that if the Lagrangean for a scalar field coupled to General Relativity only contains 

I derivatives, then it is possible to completely deparametrise the theory. This means that 

1. Physical observables, i.e. functions which Poisson commute with the spatial diffeomorphism and 
Hamiltonian constraints of General Relativity, can be easily constructed. 

2. The physical time evolution of those observables is generated by a natural physical Hamiltonian 
I which is (constrained to be) positive. 

The mechanism by which this works is due to Brown and Kuchar. In order that the physical 
Hamiltonian is close to the Hamiltonian of the standard model and the one used in cosmology, the 
required Lagrangean must be that of a Dirac - Born - Infeld type. Such matter has been indepen- 
dently introduced previously by cosmologists in the context of k - essence due to Armendariz-Picon, 
Mukhanov and Steinhardt in order to solve the cosmological coincidence (dark energy) problem. We 
arrive at it by totally unrelated physical considerations originating from quantum gravity. 

Our manifestly gauge invariant approach leads to important modifictaions of the interpretation 
and the the analytical appearance of the standard FRW equations of classical cosmology in the late 
universe. In particular, our concrete model implies that the universe should recollapse at late times 
on purely classical grounds. 
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1 Introduction 



By "the problem of time" in General Relativity (GR) one means that GR is a completely parametrised 
system. That is, there is no natural notion of time due to the diffeomorphism invariance of the theory 
and therefore the canonical Hamiltonian which generates time reparametrisations vanishes. In fact, 
instead of a Hamiltonian there arc an infinite number of spatial diffeomorphism and Hamiltonian 
constraints respectively, of which the canonical Hamiltonian is a linear combination, which generate 
infinitesimal spacetime diffeomorphisms^ . Physical observables, sometimes called Dirac observables, 
are functions on phase space which are gauge invariant, that is, they Poisson commute with all 
constraints. In particular, they do not evolve with respect to the canonical Hamiltonian. Hence 
"nothing seems to happen in quantum gravity". 

The problem of time is not only of academic interest. One of the motivations for the present 
article actually comes from cosmology and can be phrased as the following question: 

Why is it that the FRW equations describe the physical time evolution which is actu- 
ally observed for instance through red shift experiments, of physical, that is observable, 
quantities such as the scale parameter? 

The puzzle here is that these observed quantities are mathematically described by functions on 

the phase space which do not Poisson commute with the constraints! Hence they are not gauge 
invariant and therefore should not be observable in obvious contradiction to reality. Moreover, the 
time evolution described by the FRW equations is obtained from the Hamiltonian equations of mo- 
tion generated by the Hamiltonian constraint and not by an actual Hamiltonian. This is due to the 
fact that the "Hamiltonian" used to derive the FRW equations is actually constrained to vanish by 
one of the Einstein equations. The "evolution equations" generated by a constraint must therefore 
be interpreted as gauge transformations and those, by the very definition of gauge transformations, 
are also not observable, again in sharp contradiction to observation. Thus we arrive at the following 
devastating conclusion: 

Either the mathematical formalism, which has been tested experimetally so excellently 
in other gauge theories such as QED, is inappropriate or we axe missing some new 
physics. 

We will show in this article that the problem of time and the above puzzle can be solved in the 
canonical approach to GR if one manages to deparametrise the theory. By this we mean that it 
is possible to write the Hamiltonian constraints in the form C{x) = 7r{x) + H{x) where vr is the 
momentum conjugate to a scalar field 4> and where H{x) is a positive function on phase space^ which 
depends on neither or tt. In this situation it is possible to construct explicitly physical observables 
and the function H :— J d^x H{x) is the natural physical Hamiltonian which generates the time 
evolution of those observables. We will show explicitly that the scalar matter Lagrangean can be 
chosen in such a way that the physical Hamiltonian is close to the Hamiltonian of the standard model 
and the one used in cosmology and that the gauge invariant physical observables are closely related 
to the "non - observables" mentioned above. ^ 

^When the equations of motion hold. 

^As always in the canonical approach wc assume that the spacetime manifold is diffeomorphic to R, x cr where a is 
an arbitrary three manifold and x are local coordinates on a. 

^Of course it is conceivable that other than scalar matter can induce deparametrisation while the corresponding 
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The missing physics could therefore be a scalar matter component which in a precise sense is 
pure gauge. Wc therefore call it a phantom field because it is not directly observable. This is phe- 
nomenologically appealing because scalar matter has not yet been observed in nature. Its main effect 
is that it provides a notion of physical time evolution, it is a perfect physical clock. Although it is 
pure gauge, its presence has further observable consequences: The physical Hamiltonian deviates 
slightly from the usual Hamiltonian that one uses in the standard model or cosmology and therefore 
changes the dynamics slightly. The associated modified dynamics of observable quantities can be 
used in principle in order to test a given, deparametrisation generating, model experimentally. In 
fact, the modified evolution equations generated by the physical Hamiltonian rather than the Hamil- 
tonian Constraint can be recasted into FRW form, however, at least for the concrete realisation of 
deparametrisation that we consider here, now the FRW equations adopt additional terms which are 
dynamically generated. There are two types of modifications. The first type is expected: In the 
standard interpretation of the FRW equations, these can be interpreted as matter terms which at 
early times statisfy the equation of state of dust w — while at late times it becomes a cosmological 
constant w — —1. However, the energy of the scalar field is negative which requires that there be 
positive energy matter with those equations of state in order to have overall positive matter energy. 
Thus the model could be able to explain dynamically why there must be dark matter and dark 
energy. The second type of modifications are deviations from the FRW form itself. At very late 
times, where "late" depends on the parameters of the model, the FRW interpretation breaks down 
and the universe evolves drastically differently with respect to the physical Hamiltonian. In fact, our 
concrete model suggests that the universe should recollapse on purely classical grounds. Therefore, if 
we really observe evolution with respect to the physical time parameter corresponding to the physical 
Hamiltonian induced by our scalar field then the FRW equations are an approximation to the true 
evolution of the universe, which is valid at sufficiently early times of the universe only. Of course, 
the parameters of the model and its dynamical constants of motion can be tuned such that the FRW 
equations are still valid today. Let us finish this paragraph with the following provocative lesson: 

All textbooks on classical GR incorrectly describe the Ftiedmann equations as physi- 
cal evolution equations rather than what they really are, namely gauge transformation 
equations. The true evolution equations acquire possibly observable modifications to 
the gauge transformation equations vi^hose magnitude depends on the physical clock 
that one uses to deparametrise the gauge transformation equations. 

Both types of modifications just mentioned will of course not only happen in homogeneous cos- 
mology but also in full GR. Notice that we do not exclude observable scalar matter such as an 
inflaton in the Lagrangean, rather we propose that whatever scalar or other observable matter is 
present in nature, there is in addition our negative energy scalar field which is actually the reason for 
why that other matter can mathematically be related to gauge invariant, i.e. observable, quantities. 
In a sense, the mathematical formalism (gauge theory) together with the experimental evidence (e.g. 
the experimental verification of the FRW equations) inavoidably force us to conclude that there is 
something like a negative energy matter field which therefore could be called a prediction'^. 

In this paper we show that it is possible to find a whole class of scalar field Lagrangeans with the 

Hamiltonian has the properties mentioned. In particular, it would be desirable to find a scalar mode among the 
gravitational degrees of freedom leading to deparametrisation. However, this has proved to be impossible. 

^Of course, there may be other realisations of deparametrisation, different from a scalar field. However, the 
conclusion that there is a matter component of which we are unaware when we treat the FRW equations as if they 
came from a true Hamiltonian rather than the Hamiltonian constraint, remains. 
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required properties. The mechanism which leads to deparametrisation rests on a observation due to 
Brown and Kuchar made in their seminal work which enabled them to reformulate the Hamil- 
tonian constraints of GR such that they Poisson commute among each other, which is a necessary 
condition for deparametrisation as we will see. The only requirement is that the covariant scalar field 
Lagrangean depends only on the first derivatives of the scalar field. However, it may nevertheless self 

- interact due to a non - polynomial Lagrangean similar to quintessence fields [H] and more generally 
as in in k - essence models^ due to Armendariz-Picon, Mukhanov and Steinhardt All possible 
mutually Poisson commuting Hamiltonian constraints have been found in j^j, but only a subclass of 
them originate from a covariant Lagrangean which we will provide in this paper. A, possibly unique, 
two parameter family within that class leads to physical Hamiltonians which approach the standard 
model Hamiltonian when the scalar field is close to being spatially homogeneous and that of standard 
cosmology at sufficiently early times. That it is spatially homogeneous (i.e. a constant) turns out to 
be a natural requirement in order that the scalar field defines a good (i.e. synchronised everywhere 
on 0") clock. 

Curiously, as we will see, this family of scalar field Lagrangeans, to which we are driven naturally 
by physical and mathematical considerations, has been considered before by cosmologists for 
entirely different reasons. Its physical properties agree with what cosmologists call a phantom field^. 
It turns out that our family of Lagrangeans are necessarily of Dirac - Born - Infeld type with a 
constant potential. Notice again that this phantom field is not directly observable. However, we 
can, and probably must in order to have a positive energy budget, add further k - essence matter. 
Such observable k - essence matter is being discussed as a candidate for dark energy and inflation 
by cosmologists. 

We should mention that the deparametrisation technique is a special, very simple case of the 
more general "relational" approach due to Rovelli, see [7j and references therein. The mathematical 
implementation of this idea has been much improved recently |H] (see also [H]). It consists in choosing 
an infinite number of gauge fixing conditions called "clocks" and the afore mentioned physical ob- 
servables are the gauge invariant extensions, off the associated gauge cut, of non - invariant "partial 
observables" . The analytical formulae are very complicated power series in general and there are 
unsolved mathematical issues such as convergence of the series. In contrast to the deparametrisation 
case, in the more general case the associated physical observables Poisson commute only weakly with 
the constraints, that is, when the constraints hold, they are weak Dirac observables. Observables 
coming from the deparametrised theory are strong Dirac observables which is mathematically much 
more convenient. Fortunately, the much more complicated partial observable machinery is not needed 
in order to arrive at the results of the present paper. All the results that we claim in this paper will 
be proved by elementary methods, the paper is self - contained in that respect. 

We emphasise that the formalism developed in this paper is exact and non - perturbative. On 
the other hand, it is purely classical only so far. This is true for almost all the available literature 
on relational physics. In order to apply quantum theory to it, operator ordering issues have to be 
solved for the power series. This is a difficult issue in the general relational framework, however, 
under natural mathematical assumptions, we can actually solve the operator ordering problem as we 
will sketch in section |7| Yet, it may be necessary to develop a perturbative scheme just like in S - 
matrix theory. This is a good approximation as long as the (kinematical) states with which we probe 
these observables are strongly peaked, at physical time r, at the phantom field value cf) = t which 

^Basically, a k - essence field is a scalar field $ which depends non - linearly on the kinetic term g^'^ ^ ^fi^ .v ■ 
^The generally accepted rough definition of a phantom field seems to be that in a cosmological setting the first 
order term in g^'^<^^^<^ ,y of the Lagrangean comes with a coefficient which has a sign opposite to the sign in the Klein 

- Gordon Lagrangean. 
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explains why the phantom field should be close to spatially homogeneous. Hence, for a sufficiently 
short period period of physical time r, the approximation should be quite good. 

With the formalism developed in this paper, a natural platform for carrying out cosmological 
quantum field theory^ within the framework of Loop Quantum Gravity (LQG) [10] is launched. See 
[TT] for a corresponding proposal. 

The article is organised as follows: 

As this article is intended for both cosmologists and quantum geometers, in section two we state 
the results of our analysis without proofs. The proofs will be supplied in the remaining sections. 
Readers just interested in the results can therefore skip all the rest of the paper except for section 
seven. 

In section three we review the Brown - Kuchaf mechanism to generate mutually commuting 
Hamiltonian constraints. 

In section four we define the physical observables of the theory as well as the physical Hamiltonian 
originating from a general phantom field Lagrangean. 

In section five we show that physical and mathematical considerations naturally lead to a Dirac - 
Born - Infeld scalar field Lagrangean which for certain parameter range has the interpretation of what 
cosmologists call a phantom. The associated physical selection principle is that the corresponding 
physical Hamiltonian is positive and close to that of the standard model (when the metric is fiat). 

In section six we derive the consequences of the gauge invariance principle for cosmology by 
computing and interpreting the modified FRW equations. 

In section seven we conclude and outline what we plan to do with our formalism in the future, 
in particular in quantum cosmology. 

2 Summary 

The scalar field Lagrangean which leads to deparametrisation, induces a positive Hamiltonian which 
is close to that used in the standard model when the metric is fiat and which leads to physical 
equations of motion which are in agreement with the cosmological FRW equations is given by 



Here are constants of dimension'^ cm~^. We must have necessarily a > as we will see below. 
The sign of (3 is unconstrained. A natural value for (3 would be /? = a so that for small (V$)^ the 
Lagrangean becomes a;(V$)^/2 which up to the positive constant a is the massless Klein - Gordon 
field Lagrangean with the wrong sign, i.e. it is a phantom. The other natural value for /3 is = 
because (3 could always be absorbed into a cosmological term. Let us choose /? = for concreteness 
in this preliminary discussion. Lagrangeans of the form (j2.H) are being discussed in k - essence 0], 
albeit there with non - trivial potential. For our purposes, non - trivial potentials are forbidden. 

^By this we mean Quantum Gravity in the sector whose classical limit is classical cosmology. This should not be 
confused with quantum cosmology which is just a quantum mechanical toy model of the actual quantum gravitational 
field theory. 

®We assume signature (—,+,+,+) and choose units for which STrGNowton = 1- Moreover, we assume that spatial 
coordinates are dimensionless while time coordinates have dimension cm and ds^ = g^i,dX^dX^ has dimension cm^. 
We take $ to have dimension cm so that the argument of the root in (|2.1() is dimensionfrce. 



.V 
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We arrive at the model (j2.H) by a totally independent mathematical and physical reasoning, namely 
deparametrisation, hence the fact that we stumble on k - essence is rather curious. 

The canonical formulation leads to the following spatial diffeomorphism and Hamiltonian con- 
straints respectively 

= Da + T^<P,a (2.2) 



C*°* = C-y^[l + g-V,a0,fc][7r2 + a2det(g)] 

Here Da-, C respectively are the contributions to the spatial diffeomorphism constraint and the 
Hamiltonian constraint of the gravitational and non - phantom matter degrees of freedom, vr is the 
momentum conjugate to and Qab is the metric intrinsic to the spatial slices with inverse q""^ . Clearly 
a, 5, .. = 1, 2, 3 while /i, .. = 0, 1, 2, 3. From (j2.2|) we see that C is constrained to be positive. This 
will be important for what follows. If we had chosen the other sign for a then the root in ()2.2|) would 
come with the opposite sign and C would be constrained to be negative. One can also not reverse 
the sign in front of (V$)^ in ()2.1|) because this would lead to a non - definite argument of the root 
in (Q. 

The interpretation of (j2.2j) is that these are constraints, i.e. they must vanish. The canonical 
transformations on phase space that they generate are therefore not evolutions but gauge transforma- 
tions. In fact, one can show |12j that when the Einstein equations hold, the canonical transformations 
that they generate precisely coincide with spacetime diffeomorphisms. Any object which has non 
- vanishing Poisson brackets with the constraints is therefore not observable because only gauge 
invariant objects have physical meaning. The problem of time is therefore that we do not have a 
priori a Hamiltonian which generates physical time evolution of gauge invariant objects. 

The Brown - Kuchaf mechanism ^ consists in the crucial observation that 

rr\'''(l),a<l>,b = q^'DaDb =: D (2.3) 

when = 0. Thus, using ()2.3p we can solve C*°* = for tt and obtain a different Hamiltonian 
constraint 



-D- a^Q] + -D- a^Qf - a^DQ 



(x) =: ■k{x)+H{x) (2.4) 



where Q := det(g). Together with the D*°*(x) it defines the same constraint surface^ as the system 
()2.2p and is also first class. By virtue of the Brown - Kuchaf mechanism, the new Hamiltonmian 
constraints even mutually Poisson commute among each other. The arguments of the roots in ()2.4p 
are constrained to be non - negative as the derivation of that expression reveals. 

Since H no longer depends vr, 0, we have managed to deparametrise General Relativity and in 
fact the quantity 

H := ! d^xH{x) (2.5) 



^More precisely a subset of the full constraint surface. There are altogether four components of the constraint 
surface corresponding to the four possible combinations of signs in front of the two square roots involved in H2.4|l . While 
these components connect in lower dimensional submanifolds, each of them is preserved by the gauge transformations 
induced by the full Hamiltonian constraint. We therefore restrict from now on once and for all to the subset defined 
by C"°* = 0. 



6 



is a positive Hamiltonian, it is not constrained to vanish and it is gauge invariant, it Poisson commutes 
with all constraints. Next let for any real number r 



Hr := j £x [t - 0(x)] H{x) 

J a 



(2.6) 



Let / be any spatially diffeomorphism invariant quantity on phase space which does not depend on 
(j). Such functions are trivial to construct, a simple example is the volume / = / d^x^ydet{q). Then 
the series 

Ofir) :=f + {H^, /} + {H^, /}} + ^{H^, {H^, {H^, /}}} + ... (2.7) 

defines a one parameter family of gauge invariant function on phase space. Moreover, we have 

^ = WO,(r)} (2.8) 

In other words, the map r Of{T) describes the physical time evolution of gauge invariant objects 
generated by the Hamiltonian ()2.5|1 . 

The crucial additional property of H in ()2.4j) which was used in order to select the model 1)2.11] 
is that when the scalar field (j) is spatially homogeneous, which is natural in order that it defines an 
everywhere (on a) synchronised clock (f){x) = r, if the spatial diffeomorphism coinstraint holds and 
if a is sufficiently small then H{x) ~ |C(a;)| = C{x). Hence H approximates the standard model 
Hamiltonian when the spacetime is close to being flat. If we had chosen the other sign for a we 
would get H{x) ~ |C*(a;)| = —C{x). Since C{x) > for usual matter when space is flat, it would 
follow that with this sign we cannot have flat space and moreover that all matter contributions to 
the Hamiltonian come with the wrong (negative) sign. Hence the choice a > is the only suitable 
one for our purposes. 

Let us investigate this more closely: 
Since C = C^""" + C'^^^^'" + C^^-^^^^^^ = where C"^"'^'^"*^^ > is the standard matter energy 
density, we must have C = C^^'^^ + (7°s-mattor ^ g t^y^qj-q ^j^g latter contribution is from non standard 
matter such as our scalar field. The gravitational contribution C^'^^'^ is indefinite, there are positive 
and negative scalar modes contained in it and therefore it would be desirable to use a negative gravi- 
tational mode for deparametrisation. Unfortunately, such a mode does not lead to deparametrisation 
because C depends on both the corresponding vr, (p. Hence we consider (^ns-matter ^ g turns out 
that we must restrict on the portion of phase space where we have H ^ -[-[(^grav _j_ ^s-matter|^ hence a 
priori both signs in front of the square root in ()2.4|) are allowed. Hence we should have (^^is-matter ^ g 
or (^ns-matter ^ g j-espectivcly in order that H ^ (^grav _|_ ^s-matter cQ]2ies with the positive sign in 
front of (^s-matter j^gQ^use ou flat space this is the energy density of standard matter. However, as we 
will see, if (7°s-™atter ^ g ^j^gj-^ ^j^g physical evolution equations adopt modifications which lead to a 
big rip singularity in cosmological applications (the universe reaches infinite size in a finite amount 
of time). Thus, if we want to avoid this, we are naturally led to scalar matter with negative energy 
density. 

If we would have H{x) = C{x) exactly, then the physical evolution equations derived from H 
for gauge invariant observables not involving 0, vr would exactly equal the gauge transformations on 
non - gauge invariant quantities not involving 0, vr derived from the canonical "Hamiltonian" (with 
AT = 1, A^*^ = 0) 

i7™(Ar,iV)= [ d^x[N{x)Cix) +N''{x)Daix)] (2.9) 
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which ignores the phantom field. This would justify why the Hamiltonian constraint integrated 
against unit lapse is often used as a Hamiltonian. The phantom field, being pure gauge, would have 
absolutely no visible effect. However, since H{x), C{x) do not exactly coincide, there are impor- 
tant modifications, both technically and conceptually, to which we turn now in a cosmological setting. 

Namely, we will see that the FRW equations must be provided with a new and gauge invariant 
interpretation. The actual physical evolution equations generated by the Hamiltonian leads to dras- 
tic modifications in the very late universe while in the early universe (including today) they keep 
their standard form to a very good approximation, depending on the numerical value of a. In fiat, 
homogeneous and isotropic models the FRW line element takes the form ds^ = —dt'^ + a{t)'^dx°'dx^6ab 
with scale factor a and all constraints are identically satisfied due to the high symmetry of the Ansatz, 
except for a single Hamiltonian constraint 

p2 

C*°* = [-— + (A + Pm)a^] + Pphantoma^ =: C + Pphantoma^ (2.10) 

where 



7r2 



Pphantom = + =: -a v^TTx (2.11) 



is the negative phantom energy, P is the momentum conjugate to a, A is a cosmological constant 
and Pm is the energy density of all non - phantom matter. The important quantity 



X := ^ 2.12 



will be called the deviation parameter. The phantom pressure is positive 



Pphantom = "Tp? <9(aVphantom)/<9a = a (2.13) 

oa^ V 1 + a; 

leading to an equation of state and speed of sound respectively 

^ Pphantom 1 dpp]^^Yitom / dx ^2 ^2 

Pphantom 1 -|- X "^Pphantom /dx ^ • ' 



We can now do two, conceptually very different, things: 
1. 

First we follow the standard procedure in cosmology. That is, we use the constraint C*°* as if it 
was a Hamiltonian. The associated equations of motions of non - observable quantities such as the 
scale factor then lead to the usual FRW equations. From the point of view of gauge theory, the 
interpretation of those FRW equations as evolution equations of observable quantities is, however, 
completely wrong. That a{t) is not observable, that is, not gauge invariant, can be easily seen from 
the fact that da{t)/dt = {C*°*,a} 7^ 0. The correct interpretation of those equations is that they 
describe the behaviour of non - observable quantities under the gauge transformations generated by 
the Hamiltonian constraint. 
2. 

The second thing that we can do is to compute the gauge invariant functions such as Oa{T) using 
fl2.7|) with Ht- = {t — (j))H where the Hamiltonian ()2.6|) becomes 

H := VC2 - a^a^ (2.15) 
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and compute their physical evolution equations generated by (|2.15|) . 



Mathematically the two procedures are very similar to each other: In the first approach we compute 
da/dt = {C*°*,a} and express P in terms of da/dt. The first FRW equation then results by substi- 
tuting P in terms of da/dt into the constraint equation C*°*. Then we compute dP/dt = {C*°*,P} 
and insert this into d'^a/dt^ = {C*°*, a}} which results in the second FRW equation. They take 

the usual form 

3( ) = A + Pm + Pphantom (2.16) 

Oj 

d^a/dt^ _ A q 

<J — A '^[Pm + '^Pm "I" Pphantom "I" <jPphantomJ 

In the second approach we compute dOa{r) = {H,Oa{T)} = 0{H,a}{T) and can then solve Op(r) 
in terms of dOair). The first FRW equation then results by expressing C*°* in terms of physical 
observables, that is, computing Octot(r) and imposing Octot(r) = 0. That this should hold follows 
from dOc^ot{T) / dr = 0[h,c^°^}{''') = since H is an observable, hence Octot(r) = Oc^ot[(j)) = C*°* = 0. 
The second FRW equation then is obtained by computing dOp{T)/dT = 0{h,p}{'t) and using this in 

\H,{H,a}}{ 



d^Oair) / dr'^ = OsHJH,a}}{T)- This results in 



3{^^r = [A + 0,„+0,^_J(l + i) (2.17) 
3^^%^ = A(l + Vk[0p„+0,^,_J(1- V3[0,^+0,^,_J(l + i)} 
where now ^ 

and where E = H = —it is a constant of motion, namely the energy of the universe. 

Comparing ()2.17|) and ()2.18|) reveals: 
1. 

Although from the point of gauge theory it is incorrect to interpret the FRW equations (j2.16p as 
evolution equations of observable quantities, as long as x is large, the actual physical evolution equa- 
tions of observables ()2.18|) generated by the physical Hamiltonian take exactly the same form. All 
that we have to do is to make the substitution (t, a(t)) (r, Oa('r)). 

2. When x gets small, the correct equations ()2.17|) differ drastically from the incorrect equations 
(I2.16|) . Notice that what we observe in experiment is really a gauge invariant object such as Oa(r) 
and not a{t). Of course, the concrete scenario for deparametrisation that we have proposed here may 
not be realised in nature, however, we insist that whatever matter is used for deparametrisation, 
there will be corrections to the standard FRW equations. This should have observable consequences! 

Notice that we do not doubt the validity the Einstein equations fl2.16|) . They follow from the 
fundamental object C*°* which we also used in our construction. However, we stress that their 
interpretation as physical evolution equations of observables is fundamentally wrong. The domain 
of validity of the interpretation of the usual FRW equations as evolution equations is controlled by 
the deviation parameter x. It depends on the kinematical model parameter a and the dynamical 
constant of motion E. The critical value is x = 1 and is reached at scale factor Oa = \J Ejot which 
can be as large as we want for sufficiently small a. Thus we see that the mathematical formalism 
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together with our concrete model predicts that the universe evolves differently at late times, that is, 
at large scale factor. We expect similar modifications in other applications of GR such as black hole 
physics and it is an interesting speculation that the corresponding gauge invariant interpretation of 
Einstein's equations could predict large scale deviations from Newton's law which then could be in 
agreement with the measured rotation curves of galaxies. Notice that all of this is a purely classical 
effect, there is no quantum gravity involved in this although our motivation, deparametrisation, 
certainly comes from quantum gravity. 

The fact that the phantom makes a negative contribution to the energy budget may be disap- 
pointing for supporters of k - essence where > is usually required. However, Pphantom < is of 
no concern as long as the remaining matter makes an overall positive contribution^". Actually, since 
the gravitational contribution to the Hamiltonian constraint is negative definite in cosmology, in fact 
due to C*°* = we must have pm > |pphantom|- Notice by the equation of state the phantom behaves 
like dust at small scales Pphantom — > ~E/0^ and as a negative cosmological constant Pphantom — —ol 
at large scales. This can be easily compensated by additional positive energy k - essence matter 
or simply by ordinary (dark?) matter plus an additional cosmological constant term A — a > 0. 
In a sense, if we want to explain the observational fact that the FRW equations describe the uni- 
verse while their mathematical derivation violates the principles of gauge theory, then something 
like a phantom is needed for deparametrisation and in turn it requires something like k - essence 
for reasons of total positive matter energy budget. From this point of view, both a phantom and k 
- essence are a prediction of the mathematical formalism (gauge theory) together with observation 
(FRW cosmology). 

We will see furthermore that in order that the universe does not reach infinite size in finite r 
time, it must in fact recoUapse which can be achieved by a suitable choice of the parameters. Then 
the picture of a periodic universe arises if one can establish that Quantum Gravity effects avoid big 
bang and big crunch singularities. This would imply that the universe evolves through the "would 
be" singularities in an infinite number of cycles. Notice that recently ^3] a simple cosmological toy 
model has been rigorously quantised by the methods of LQG using precisely the gauge invariant 
programme suggested in jTT] and for which possible classical foundations have been layed out in the 
present paper for the full theory. In that model, the singularity is indeed quantum mechanically 
avoided which is a promising hint that in full LQG the singularity is avoided as well. 

We close this section with some final remarks: 
1. 

From the point of view of a cosmologist nothing would be more natural than to use the scale factor 
itself as a clock: It is a monotonic function of the unphysical time parameter t (until possible rec- 
oUapse). Why did we not do that immediately (we can do it indirectly, see below)? There are two 
reasons. First of all, we wanted to provide a universal framework, i.e. to provide a physical notion 
of time in all possible situations and not only in homogeneous ones. However, in inhomogeneous 
situations, the notion of a scale factor is void. As a substitute one could consider the volume of 
(subsets of) a (the total volume is infinite for non compact topology of cr). However, this does not 
work for the same second reason for which also the scale factor itself is inappropriate in the homoge- 
neous situation: In order to achieve deparametrisation and to obtain a physical Hamiltonian with all 
the required properties, the clock variable (s) must be cyclic in the Hamiltonian constraint (at least 
weakly). This condition is violated for the scale factor and its inhomogeneous relatives due to the 
universal coupling of matter to gravity. 

^°The usual energy conditions on the energy momentum tensor do not make any restrictions on individual matter 
species but only on the overall matter content of nature. 
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2. 

Of course, after one has deparametrised the system and one is only dealing with physical quantities 
such as the physical scale factor Oair) (or one of its inhomogeneous relatives) one can use it as a 
physical clock in place of r itself which is maybe better geared to what one does in reality. One can 
then express time dependence of other observables Of{T) in terms of Oa by solving Oa(r) for r. In 
other words, while we cannot use the scale factor to deparametrise the system, we can still use it as 
physical clock after deparametrisation. The time evolution in terms of the physical scale factor will 
then be generated by a more complicated physical Hamiltonian. 
3. 

One could think that what cosmologists usually do in order to describe measurable quantities math- 
ematically is actually precisely correct, that is "relational". For instance the redshift factor 



,(f n - ~ ""^^'^ 

Z[ti,t2 .— 7— 

iU2 a{ti) 



(2.19) 



is the ratio between the emission frequency ui of a spectral line (known from a table top experiment 
on Earth) and the absorption frequency UJ2 observed on Earth is certainly measurable. Formula 
()2.19|) relates this observable quantity to the ratio of the scale factors at unphysical emission time ti 
and absorption time t2 respectively. We will now show that ()2.19|) is in fact incorrect: 
The reason is that the quantities a{t) are not observable. In order to see what is going on, we have 
to go through the derivation of the redshift formula. Consider a star at comoving distance r from 
Earth. For light the geodesic is null and due to ds"^ = —dt^ + a{t)^dx'^dx'^6ab we get as an equation 
of motion a{t)r{t) = 1. Formula ()2.19j) then results from the fact that the beginning and the end 
of the wave travel the same comoving distance r = J^^ dt/a{t) = J^^^^^ dt / a{t) with uoj = 2tx /Tj. 
This is certainly mathematically correct, however, the quantities a{tj) are not observable. In order 
to express z in terms of observale quantities Oa{T) we express the line element in terms of r (see 



ds^ = -dr^n + -) + OaiT^dx^dx^Sab 

X 



(2.20) 



Notice that r is no gauge parameter but a physical observable associated with the physical Hamilto- 
nian, hence the factor 1 + l/x cannot be transformed away by a diffeomorphism r 1— > ip{T) without 
changing the Hamiltonian. We now obtain the null geodesic equation of motion Oa{T)dOr{T) / dr = 
V^l + l/x. The same argument now leads to the modified redshift factor relation 



Z[Ti,T2 



Ll)2 



Oa{r2) 

Gain) \ 



1 + 



x(r) 



x{t2) 



(2.21) 



and now all displayed quantities are observable. Hence we see that as long as x is large, ()2.21|) and 
fl2.19|) agree in the following sense: What one incorrectly does in cosmology is to identify the unob- 
servable gauge pair (t,a(t)) with the pbservable physical pair (r, Oa(r)). With this interpretation, 
the wrong relation ()2.19|1 is a good approximation to the correct relation 1)2.211) as long as x is large. 
However, there are large deviations especially in the late universe and of course the modification 
fl2.21|) may have an observable effect on the interpretation of supernovae type la observations (stan- 
dard candles) which provide evidence for recent accelarated expansion of the universe. 



We now proceed to the mathematical and physical details. 
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3 Review of the Brown — 



Kuchaf Mechanism 



3.1 Covariant, Minimally Coupled, Potential Free Scalar Fields 

In order to prepare for the explanation of the Brown - Kuchaf mechanism we review here the canon- 
ical formulation of general scalar field Lagrangeans of a special class. 

We consider a general, covariant scalar field Lagrangean minimally coupled to the metric with action 

Vantom = / d^X ^ | dct ((7) (7^'^$,^$,,/2) (3.1) 

where L is an arbitrary function of the variable indicated. It will be crucial for the Brown - Kuchaf 
mechanism to work that the scalar field $ only appears with derivatives, i.e. there is no non - trivial 
potential term. Obviously, ()3.1|) is invariant under Diff(M). 

As usual we perform a 3+1 split of the action jT2| and assume that M is diffeomorphic to 
R X cr where a is a three - manifold of arbitrary topology. Hence, there is a foliation t i— > = 
Yt{a) of M by spacelike hypersurfaces which are the images of a under a one parameter family 
of embeddings t Yf. This way we obtain a diffemorphism IH x cr — > M; {t,x) X := Yt{x). 
We consider the fohation vector field T{X) := [dYt{x) / dt]Yt{x)=x which can be split as T{X) = 
[N(t,x)n{X) + N°-(t,x)dYt{x)/dx°']Yt{x)=x- Here x", a = 1,2,3 are local coordinates of a while 
X^^, fi = 0, 1, 2, 3 are local coordinates of M. The vector field n is everywhere normal to the 
foliation, that is, g^yU^Y^^ = g^yn^rf + 1 = 0. The functions A^, N°- respectively are known as lapse 
and shift functions. 

We now pull back ()3.1|) by the diffeomorphism Y and express everything in terms of 
N{t,x), N"-{t,x), qabit,x) = {Yt*g)abix) and 0(t,a;) = (F/$)(a;). It is not difficult to check that in 
the embedding coordinates the components of the metric tensor read gu = —N"^ + N"'N^qab, gta = 
qabN\ gab = gab and for the inverse c/" = -l/N\ c/*" = A^VA^^ = g"'' - N'^N^/N'^ where 
(f'^gcb = ^b- It follows that (the lapse is assumed to be everywhere non - negative) 

V|det((7)| = N^det{q), I := -g^'^^^^^u = (V„0)' - q"' <p ,a(t> ,b (3.2) 

where n = (T - N'^Yt^a) /N so that V„0 = (0 - A^'^^.a) /N. 

We are now in position to perform the Legendre transform. We have 



7r(t, x) := <5S'phantom/50(t, x) = ^d^[V„0]L'(//2) (3.3) 
where the prime denotes the derivative with repspect to 1/2. From ()3.3|) we infer 

K := [-J^f = [L'{imf{I + Vl V:= g-^V^,. (3.4) 
Vdet(g) 

We assume that L is such that ()3.4|) can be solved uniquely for / = J{K,V). Then ()3.3p can be 
solved for V„0 

V„0 = p/L\J/2), p := 7r/v/det(g) (3.5) 
We can now complete the Legendre transform 

Vantom = [ dt [ d'x{7r^ - [AT'^TT^,, + AT v^dct (g) {^V^' ( J/2) - L{J/2)}] (3.6) 
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From (j3.6|) we read off the contributions of the scalar field to the spatial diffeomorphism and Hamil- 
tonian constraint respectively 



^phantom _ ^^^^ ^37^ 

3.2 The Brown — Kuchar Mechanism 

Let us denote by Da, C respectively the contribution to the spatial diffeomorphism and Hamiltonian 
constraint respectively of the gravitational field and all other matter fields (say of the standard 
model or one of its supersymmetric extensions). Then the spatial diffeomorphism constraint is given 
by = Da + and the Hamiltonian constraint by C'°^ = C + Cphantom^ rj.^^ simple, 

but crucial observation due to Brown and Kuchar is that we may use the spatial diffeomorphism 
constraint in order to remove the dependence of on altogether, thus making it a function of 
p, the gravitational field and all the other matter fields only. Namely, we have, when = 

V = q = = = =■ d/K (3.8) 

This is the Brown — Kuchar Mechanism: The field 0, which appears only in the combination V 
within the Hamiltonain constraint, has been eleminated. This would not work if the Lagrangean also 
would depend on explicitly (potential term), not only through the combination / which involves 
only derivatives of 0. 

Consider now the function J{K, d) := J{K, V = d/K). Then, the Hamiltonian constraint can be 
equivalently described by the function 

K 



= v/d^[c + l—j— - L{J/2)] (3.9) 
L {J / Z) 



where C = a/ det(g)c. Since the constraints form a first class system, also the new constraints do 
and they define the same constraint surface. 

Notice that ()3.9|) depends on p only through K. We will now assume that we may solve ()3.9p for 
K algebraically (possibly with several branches) 

K = G{c,d) (3.10) 

Notice that by construction G is (constrained to be) non - negative. We may therefore write the 
Hamiltonian constraint in the still equivalent form 



C"*°* = 71+ ^det{q)^/Gic, d) (3.11) 

The other sign is also possible but the above choice leads to a positive physical Hamiltonian close to 
that of the standard model when the metric is flat (as mentioned in section ()2.2j) we will restrict to 
the subset of the constraint surface defined by 1)3.111) in what follows). 

What is remarkable about the functions [G{c, d)] [x) is that they mutually Poisson commute among 
each other. The formal proof is as follows: The constraints ()3.11|) form a first class system, hence 
their mutual Poisson brackets is a linear combination of the constraints C"*°* = tt + a/'G(c7^, -0^°* 
with structure functions. However, since the vr Poisson commute among themselves as well as with 
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H = ^ det {q)G because G does not depend on 0, it follows that {C'^°\x), C"*°*(|/)} = {H{x), H{y)} 
does not depend on tt, (j) any more. Thus, the linear combination of constraints with structure 
functions, which are non - singular on the constraint surface, must conspire in such a way that the 
dependence of the bracket on vr, drops out completely. Suppose then that 

{C"°\N), C"°\N')} = I d'x[h,N'{x)C"°\x) + f%^^,{x)Dl"\x)] (3.12) 

where N,N' are test functions and C{N) = J d^xN{x)C{x). Since fl3.12j) does not depend on vr, 
we may choose tt such that (jH.llj) vanishes. Then only the second term in ()3.12|1 survives and must 
no longer depend on (p. It follows that 



d'x[r^^^,ix)D':\x)l_^^^^^ = j d'x[r^^^,{x)l_^^^^iD, - ^detiq)G<P,a){x) (3.13) 
We can now expand = (/^ A^')^-_^ det(g)G powers of 0,a, that is 

oo 

^" = ^?o + E<'"''"'^.^i-'^.^" (3.14) 

n=l 

where the coefficients are supposed to be independent of vr, 0. The resulting recursion relation is 
then given by (tt = — A/det(g)G being understood) 

I^a^n'"-'" + ^^n-l" = (3.15) 

and can be solved for instance by 

gn'-'" = i-T^r-'h'h'^.h'-^g'," (3.16) 

where /i" is an arbitrary function such that h'^Da 7^ and Qq is also arbitrary. However, then 

is singular on the constraint surface and in fact g°'{Da + 7r0.a) = (^g-Da is not a linear combination of 
constraints. 

We do not need to rely on such a formal argument: The rigorous proof is by actually com- 
puting the Poisson bracket. In [S] we find the necessary and sufficient condition for expressions of 
the form H = a/ det(g)G(c, d) to be mutually Poisson commuting: Consider functions of the form 
Hn{Q,c,d) := Q''/'^K{c,d) where Q = det(g). Then, using the well known Poisson algebra gener- 
ated by the Hamiltonian and spatial diffeomorphism constraints ^2] one can compute the Poisson 
brackets between the smeared functions Hn{N) := j d^xN{x)Hn{x) and ask for the condition on hn 
such that {Hn{N), Hn{N')} = 0. This leads to the following first order partial differential equation 

2^"9J = ^[^] -it^] ^'-'^^ 

Dividing this equation by h"^ we get a PDE for ln(/i„) which one can solve by the method of separation 
of variables. The general solution, also called complete integral in the theory of first order PDE's, is 
given by the two parameter family 

Tl S S — Tl 

\n\hn{c, D; a,b)] = b + - ln{d) + 2eac + -[2s + n ln( )] (3.19) 

4 4 s + n 
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where s = ^/nP~^^T6a?d and e,6 = ±1. The so called general integral, which depends on an arbitrary 
function g, is obtained by solving the equation 

d\n[hn{c, d;a,b = g{a))]/da = g'{a) + 2ec + 6^ = (3.20) 

for a in terms of c, d and to reinsert the solution a = A{c, d) into that is, /i^(c, d) := /?.„(c, d;a = 
A(c, rf),6= (7(A(c, c/))). 

Of course, in our case we are interested in the case n = 1 and all we have to do is to check whether 
hi{c,d) = \/G{c, d) solves ()3.18p . The variety of solutions h\ is certainly bigger than those that 
come from a covariant Lagrangean, that is, from a given function L. 

The reason for why we mention the variety h\ is that we will be interested in solutions h\ to ()3.18p 
with special properties and instead of guessing a suitable Lagrangean it may be more constructive to 
select a candidate in the family h\ and to ask from which Lagrangean, if any, it derives. In order to 
answer this question, suppose we have selected a solution (c, d) = G{c, d) of ()3.18p . Then we solve 
fl3.10|) K = G{c, d) algebraically for c = c{K, V) using again the key identity d = V/K. From this 
we infer -c{K, V) = K/L'{J/2) - L{J/2) where I = J{K, V) is the solution of K = L'{I/2f{I + V^). 
Hence we get 

c{L'{I/2f{I + V), V) = L'{I/2){I + V)- 1(1/2) (3.21) 

In order that this be an ordinary first order differential equation for L the explicit dependence of 
()3.21|) on V must cancel. In order to achieve this, one performs algebraic manipulations to 1)3.211) and 
obtains an equation which is a polynomial in V. The coefficients of all powers of V in that polynomial 
then have to vanish. This leads in general to more than one ordinary differential equation for L which 
are contradictory if L does not exist and which are equivalent if L exists. This is the necessary and 
sufficient condition for a given solution ^G{c, d) of ()TTR|l to come from a covariant Lagrangean of 
the form L(J/2). 

4 Deparametrisation of General Relativity 

The two ingredients that we need here are the following properties of the functions H{x) = ^jQG{c, d){x) 
derived in the previous section: 

1. they are mutually Poisson commuting. 

2. they do not depend on tt, 0. 

3. they are densities of weight n = 1. 

This is all we need in order to show that the following object 

Ofir) := K/(/)]m=.-<^, «m(/) = n^^^^^^ ^^'^^ 

n=0 

where H{M) = d^xM{x)H{x) for some smearing function M, is a one parameter family of strong 
Dirac observables, i.e. it Poisson commutes with both the spatial diffeomorphism constraint and the 
Hamiltonian constraint. Due to the Poisson commutativity of (f){x),H{y), 1)4.11) can also be written 
as 

oo „ 

Ofir) = ar-M = E '■= / ^'^[^ - <P{x)]H{x) (4.2) 
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The multiple Poisson bracket appearing in (j4.1|) . (j4.2j) is defined iteratively by {i^T-, /}(o) = /, 

{Ht, /}(n+l) = {Hr, {Hr, /}(n)}- 

The requirement on / is that 

A. it is already spatially diffeomorphism invariant^^ and that 

B. it does not depend on vr, 0. 

For the expert the invariance of Of{T) under the gauge motions of GR should be already obvious: 
First of all, by inspection, H.^ is spatially diffeomorphism invariant because [(r — (j))H]{x) is a scalar 
density of weight one. Hence Of{T) is spatially diffeomorphism invariant. Next, let us consider the 
Poisson automorphism f^M^f) '■= exp(X(7tot(M)) ■ / where Xp denotes the Hamiltonian vector field 
of some phase space function F and C*°*(M) = J^d^xM{x)C^"^{x) for some test function M. Then 
we have 

MOfir)) = PM{a^-Af)) (4.3) 

= a/3„(r-9i)(/5Af(/)) 

= ar-^-M{PM{f)) 

= ar-0-M(«Af(/)) 

= ar^M = Of{r) (4.4) 

where in the second step we used exp(X) exp(y) = exp([exp(X)y exp(— F)] exp(y) for Hamiltonian 
vector fields X,Y, in the third we evaluated (3M{Hr) = Hr — H{M), in the fourth we noticed that 
/3m = C(M on functions which do not depend on and in the last we noticed that Hr, H{M) Poisson 
commute so that a^-^-M = ctr-^ ° c^m • 

The next two subsections can therefore be skipped by the expert, for the non - expert we will 
supply more details about symplectic geometry and will carry out the calculations in a less elegant 
but more pedestrian way which has the advantage of explicitly displaying the mechanism due to 
which all of this works. 



4.1 Invariance under Spatial Diffeomorphisms 

In order to establish spatial diffeomorphism invariance, let us write ()4.H) more explicitly as 

Of{r) = f + Y.-\ I ^'^1- / - - 0(x„)){i7(xi), ..{i/(x„,), /}..} (4.5) 

where we have made use of the fact that 0(x) Poisson commutes with H{y), f by assumption. 

We will establish a result which is stronger than mere Poisson commutativity of Of{T) with 
D^°^{u) = J d^xu°'{x)D^°^{x) for arbitrary vector fields u on a: Let s i— > be the one parameter 
family of spatial diffeomorphisms defined by the integral curves c"(s) which are the unique solutions 
of the system of ordinary differential equations c^(0) = x, c"(s) = m(c"(s)), that is: ip'^{x) := c"(s). 
We define for arbitrary functions / on phase space 

n=0 

^^Such functions are easy to construct, any integral of a scalar density of weight one constructed from the canonical 
variables is spatially diffeomorphism invariant. 

^^For completeness, let us mention |H| El that if / does depend on either tt or then (|4.1() must be generalised to 
Otir) :— [(3M{f)]M=T-(t>- f independent of TT,(j) we have /3Af(/) = ctMif), hence we arrive at (|4.1|l . 
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Obviously {^)s^Qa^u{f) = It follows that s i— >■ a^u is a one parameter family of 

Poisson automorphisms^^ with the spatial diffeomorphism constraint as generator and one easily 
checks that a^u(T — (j){x)) = r — (f){(p^{x)) and a^u(H{x)) = \ det{d(p^{x) /dx)\H{(p^{x)). Hence, 
the automorphisms generate transformations on phase space which are in agreement with the fact 
that under a spatial diffeomorphism (f G Diff(cr) the function r — (f){x) i— > r — (j){ip{x)) transforms 
as a scalar because t is a constant. This would not hold if r would be a non - trivial phase space 
independent function. This will be important for our discussion below where we derive the relational 
origin of 0/(r). On the other hand H{x) i— > | det{dip / dx)\H {ip{x)) transforms as a scalar density of 
weight one. 

Since / Poisson commutes with by assumption we have Q;,^«(/) = /. Since a^pu is a 

Poisson automorphism it follows that 

a^n{[T - 0(xi)]..[r - /}..}) (4.7) 

= a^^ir - 0(a;i))..a^|(r - 0(a;„)){a^«(iy(xi)), 

= I det(9¥.:(xi)/9xi)|..| det(a(^:(x„)/9x„)|..[r - 0(y.^(xi))]..[r - 0(</^:(x„))] x 

x{i/((^:(xO),..{i/(¥^r(x.)),/}..} 

Here we have used that the Jacobean | det(9v9/9x)| commutes with Poisson brackets. Now invariance 
of ()4.5|) trivially follows since the Jacobean allows us to change variables under the integrals. 

This holds for the infinitesimal diffeomorphisms which are generated by and which allow us 
to explore the component of the identity of Diff((T). However, the invariance certainly extends to the 
full diffeomorphism group if 0, H transform as scalars of density weight zero and one respectively. 



4.2 Invariance under the Hamiltonian Constraint 

Consider the smeared Hamiltonian constraint C"*°*(M) = (i^xM(x)C"*°*(x) where C"*°*(x) = 7r(x) + 
li{x) and M is an arbitrary test function. Using the explicit expression for Ojir) in the form ()4.5|) 
we notice that there are contributions to C'^'^^{M).,0 fir)^ from {C"*°*(M), = M{x) and 

{C"*°*(M), /}..}} = I SxM{x){R{x), {Hix,), /}..}} (4.8) 

where we have used the fact that /, H{x) and therefore {if (xi), ..{if (a;„), /}..} do not depend on cj). 
We compute 

{C""\M), Ofir)} = {HiM), /} + 1 /■ dW. [ d'xn x (4.9) 

n 

[- - Hx,)]..[r - 0(xfc_i)]M(x,)[r - 0(xfc+i]..[r - 0(x„)]{ii(xi), ..{ii(x„), /}..} 



X ^ 

k=l 

+ [ d'xoM{xo)[r-<P{xi)]..[T-<P{xn)]{H{xo),..{H{xn),f}..}] 



The crucial Poisson commutativity of the H{xk) implies that the function 

{xi,..,Xn) ^ {H{xi), ..{H{xn), f}--} is Completely symmetric in its arguments^'^. Thus, relabelling 



i^That is, {«(/), a(g)} = a{{f,g}) and a{f + g) = a{f) + a{g),a{fg) = a{f)a{g), a{f) = a{J) for arbitrary, 
possibly complex valued, phase space functions f,g. 

i^Proof: We have {iJ(Mi), ..{iJ(Af„), /}..} = X^(Mi) • ^//(a/.) • ■• • ^H(Af„) • / where H{M) = J d^xM{x)H{x) 
with arbitrary test functions A/j, and X^f(A/) is the Hamiltonian vector field of H{M) which acts on fimctions by 
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x'q := Xki x\ := Xi, .., x'i^_^ := Xk~i, x'f. := Xk+i, := Xn we can write (|4.9|) in the form 

{C'''^\M),OfiT)} = {H{M)J} 

- Yl ( \a| / ^^^O-- / d^Xn-lM{xQ)[T - <P{xi)]..[t - (p{Xn-l)]{H {xq) , ..{H{Xn-l), /}••} 
n=l ■^'^ 

+ / ^'^O- / M{xo)[T-^{x,)]..[T-4>{Xn)]{H{xo),..{H{Xn),f}..}] 

n=l •^'^ •^'^ 

= (4.10) 
which finishes the proof. 

4.3 Physical Hamiltonian 

We claim that 

H := [ d^xH{x) (4.11) 

is a physical Hamiltonian, that is, a physical observable which generates the time evolution r 0/(r) 
of all physical observables. 

That H is spatially diffeomorphism invariant is trivial because it is the integral of a scalar density 
of weight one. That it Poisson commutes with all the C"*°*(M) is also trivial because the li[x) 
mutually Poisson commute among each other. Hence if is a physical observable. 

To see that H generates the r evolution we compute 

— 0/(r) = ^ — ^ / d^xi.. / c?^x„[r - 0(xi)]..[r - 0(xA:-i)][r - 0(xfc+i)]..[r - 0(x„)] X 

n=l ^' fc=i J'^ J'^ 

x{ii(a;i),..{if(x„),/}..} 

= J ^^^0- / " (p{Xi)]..[T - (f){Xn-l)]{H{Xo), {H{X - 1), ..{H{Xn), /}••} 



T - 0(xi)]..[r - <P{x^)]{H{x - 1), ..{Hixn), /}..}} 



= {H.Oj{r)} (4.12) 

where we have made use of the same manipulations as in the previous section and in the last step 
we used again that {H{x), (f>{y)} = 0. From ()4.12p it follows in particular that 

oo „ 

Of{T) = aAOf), Of := Of{0), ar{F) := ^ ^{H, (4.13) 

n=0 

which can be checked explicitly by expanding Of{T) in powers of r. In practice one therefore computes 
Ofir) by solving df{T)/dT = {if, /(r)} and then chooses the "constant" of integration to be such 
that /(r = <P) = f. 

Xh(M) ■/ = {H{M), /}. The relation [Xh(m)- Xh{m')\ = X[h{m).h(M')} holds between the Lie bracket of Hamitonian 
vector fields and the associated Poisson brackets. The claim now follows since {H{M), H{M')} = 0. 
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Finally, let us note the the identity 



{Oj{T),OfiT)} = 0{fj,}{T) 



(4.14) 



which immediately follows from the fact that Of{T) is the Hamiltonian flow of / generated by H^-. 
Hence the physical observables satisfy a simple Poisson algebra if the /, /' do. Mathematically 
speaking, / i— > 0/(r) is a one - parameter family of Poisson homomorphisms. 

4.4 Relational Origin of the Formalism 

The following section unveils the relational origin of our formalism and can be skipped by the reader 
merely interested in the physical application of the phantom field. On the other hand, one learns in 
this section why r = const, rather than r(x) is natural and why if is a natural Hamiltonian. We 
also finish this section with some cautionary remarks which list some assumptions of the formalism 
which were not yet explicitly mentioned. 

Given a system of first class constraints Cj possibly with structure functions, suppose we find func- 
tions Tf such that the matrix A with entries Au := {C/, Tj} is invertible and let B = A~^. Consider 
the functions C'j := J2j BijCj, fix real values tj in the range of Tj and let 



where Xj denotes the Hamiltonian vector field of Cj. One can show that the Xj are weakly commut- 
ing, hence the sequence of the application of the Xj in ()4.15|1 is irrelevant on the constraint surface 
jHj. One can show that ()4.15|) is a weak Dirac observable, i.e. it Poisson commutes with all the Cj on 
the constraint surface. One can also show [H] that the evolution in tj has a Hamiltonian generator 
Hi{{tj}) which is defined via {Hi{{tj}), FJj,} ■= dFJrp/drj for those / which have vanishing Poisson 
brackets with the T/ and their conjugate momenta. However, these Hamiltonians are not granted to 
be either positive or independent of the tj ||9j. The physical meaning of FJrp is that it is the value 
of / in the gauge^^ when Tj assumes the value tj. See ||8j for more details. 

In General Relativity the label set of the /'s takes countably infinite cardinality and there are 
open issues with the convergence of ()4.15|) . In particular, the fact that ()4.15p is only a weak Dirac 
observable and Hj{{tj}) a weak Hamiltonian is mathematically rather inconvenient. Moreover, the 
inversion of the matrix A which is required in order to compute FJj, at least order by order is 
practically difficult for general choices of the Tj which is why it is important to supply physical 
input towards choosing the "clocks" Tj. Notice also that in General Relativity the Cj will involve 
both the spatial diffeomorphism and the Hamiltonian constraint even if / is spatially diffeomorphism 
invariant because the Hamiltonian constraints do not close among themselves, they are proportional 
to a diffeomorphism constraint. This can be circumvented when the Tj themselves are also spatially 
diffeomorphism invariant 8j, however, it is difficult to choose an algebraically independent set of 
such functions which also satisfy the requirement on A and which can be considered as canonical 
configuration variables. 

These comments reveal that ()4.15|) is practically difficult to handle unless one manages to simplify 
it drastically, in particular the matrix A should be simple. We claim that this is precisely what we 
managed to do in this paper: Let bj be an orthonormal basis of L2{cr, d?x) such that 6/ also has finite 
Li(R, d^x) norm^^. Let x ^ t{x) be an arbitrary function and define Tj :=< bj, r >, Tj :=< bj, > 

^^Notice that Tj = tj can be considered as a gauge fixing condition. 
^^For a — R^ these could be Hermite functions. 




(4.15) 
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, TT/ =< bj, TT >, C'j :=< bj, C"*°* > where we assume that r, 0, vr, C"*°* have at least finite sup norm 
on a. Then we find for our Ht that if^ = ^/(t/ — Ti){Ci — tt/). Notice that A/j = {C'j.Tj} = Sjj 
is the unit matrix so that C/ = C'^. Therefore, our Of{T) coincides with FJrp up to the fact that 
we use r = const, rather than arbitrary r(x). That FJrp Poisson commutes (even strongly) with the 
Cj is now a consequence of the fact that our Cj do close among themselves, namely they form an 
Abelean subalgebra of the constraint algebra. 

What is unclear however from the general relational framework is that this FJj, should Poisson 
commute with the spatial diffeomorphism constraint because formula ()4.1|) does not involve the spa- 
tial diffeomorphism constraint at all and the Tj are not at all spatially diffeomorphism invariant. 
Since FJj, coincides with Of{T) with Hr replaced by J^d^x{T{x) — (f){x))H{x) we see that this is 
automatically the case if and only if t = const.. This explains why, for gauge invariance reasons, 
r = const, is the only reasonable choice for a function and why^^ tj = t J d^xbj is only a one pa- 
rameter family of time evolutions rather than a "bubble time evolution". The generator of this r 
evolution is then 'Y^^ < l,bj >< bj,H >= H{1) = H which now has the advantage not to have an 
explicit T dependence and which is positive by our construction. 

Thus we see that, by mere coincidence, the Brown - Kuchar mechanism helps to drastically sim- 
plify the relational framework. Rather than computing the infinite number of series ()4.15|) and 
inverting complicated infinite dimensional matrices there is only one series ()4.H) and no matrix to 
invert. There is a distinguished notion of time generated by an invariant, positive and time indepen- 
dent physical Hamiltonian and the observables that we compute are strong observables, they have 
vanishing Poisson brackets with all constraints everywhere in the phase space, not only the constraint 
surface. Since there is only one series, convergence issues and, in quantum theory, operator ordering 
issues are much easier to settle. 

Remarks: 

1. There is an issue which we have not considered so far: The functions H{M),C'^°^{M),7c{M) 
may not converge for M = t = const, if a is not bounded and/or may not be functionally 
differentiable if a has a boundary. In fact, from the case of a massless Klein Gordon field which 
corresponds to the choice L{I/2) = 1/2 and in case that the geometry is asymptotically fiat 
one would assume that vr and therefore necessarily H decay only as with respect to an 
asymptotic radial coordinate and thus Hr would blow up. We will assume that this problem is 
absent by a judicious choice of phantom model with corresponding fall - off conditions for both 
vr and H. For instance we could use compact cr without boundary and these issues would be 
absent. This is not sufficient for asymptically flat boundary conditions or in = 0, 1 cosmology. 
In the model we derive in the next section this problem will be avoided automatically because 
the Hamiltonian turns out to have compact support, at least when the scalar fleld is close to 
homogeneous and the spatial diffeomorphism constraint holds. 

2. One may worry that due to square roots which enter into the construction of C"*°*(M) its 
Hamiltonian vector fleld is singular or zero on the constraint surface. Indeed, that would be 
the case if one would drop the phantom fleld altogether and consider H{x) as an alternative 
choice of the Hamiltonian constraint for GR and the remaining matter flelds which together 
with the spatial diffeomorphism constraints Da{x) actually form a Lie algebra without structure 
functions j5j. However, we do not drop the phantom fleld and due to the term ti{x) = C*°*(x) — 

^^This explains why we need < oo. 
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H{x) the Hamiltonian vector field does not vanish on the constraint surface and it is not singular 
because H{x) is not constrained to vanish, rather 7r(x) = —H{x). 

3. The observables that we have constructed here are strong Dirac observables with respect to 
the new Hamiltonian constraints C"*°*(a;) but since we have used the constraints in order to 
transform between this constraint and the original one, the observables will be only weak Dirac 
observables with respect to the original Hamiltonian constraints C*°*(x). 

4. We use the phantom field values (f){x) as physical clocks. Under the gauge transformations 
generated by the new Hamiltonian constraint C"*°* it transforms as S(j){x) = {C"*°*(M), (j){x)} = 
M{x). Since the lapse function M is required to be everywhere positive, it follows that under 
this unphysical time evolution (f){x) evolves strictly monotonously for all x. Therefore (f){x) is 
classically a perfect (i.e globally on phase space) clock for all x. 



5 Physical Scalar Field Models: 

Dirac — Born — Infeld Phantom k — Essence Lagrangeans 

In the previous sections we have outlined why particular scalar field models in general allow to 
deparametrise GR. What is left to do is to exhibit a (set of) model(s) which leads to a Hamiltonian 
that reduces to that of the standard model when the geometry is fiat. This is the task of the present 
section. The analysis displayed here is far from complete and we will satisfy ourselves by finding one 
suitable model. More general or improved models are left for future research. 



5.1 Selection of the Model 

The class of models we will look at is already constrained by the necessity to be able to solve ()3.4|) 

K=[L'(//2)f (J + r) (5.1) 

algebraically for / = J{K,V). Hence (jS.lj) should lead to an algebraic equation at most of fourth 
order in /. We will restrict attention to models which only lead to quadratic equations in order 
to avoid the algebraic complications associated with Cardano's and Ferrari's formulas for cubic and 
quartic equations resepctively. This restricts us to functions L that satisfy 

™r - (5.2) 

where a, b, 6, ^, ( are real constants such that the square roots that enter the integral L of ()5.2j) . when 
evaluated at the solution / = J, of ()5.1|1 have positive arguments. 

The integral of ()5.2j) can be carried out for all values of a,b,S,^,( but in general involves com- 
plicated inverse trigonometric functions and logarithms. This would be bad because we also need 
to solve the Hamiltonian constraint for K later on and such functions would lead to transcendental 
equations which cannot be solved algebraically. In order to avoid transcendental equations we have 
to specialise (j5.2|) to one of the following cases: 

ii. [L'{I/2)]^ = a" 

iii. [L'{I/2)f={^bf[a + bI/2] (5.4) 
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which are readily integrated to 



11. 



111. 



L(J/2) = -(/? + e^a + hi/ 2) (5.5) 
L(//2) = -p + aI/2 

L{I/2) = -{(3 + e^a + hi 12) (5.6) 

where a, 6, /3 are real constants on which we will impose some restrictions in what follows and e = ±1. 
Remarkably, the Lagrangean i. in ()5.5p is precisely of the form of a Dirac - Born - Infeld type of 
phantom field with constant potential if we set e = —1, see e.g. the first reference in jSj and references 
therein. 

Let us first discuss case i. which turns out to be the right choice. Equation ()5.1|) now leads to 



K 



4(a + hi 12) 



It follows that 



a + 6J/2 



J{K,V) 



a - hV/2 



h^ - 2Kh 



(5.7) 



(5.^ 



1 - 2K/h 

which should be positive independently of the range of both K, V. Since K, V are manifestly positive, 
in order that the square root in ()5.5|1 is well defined (real valued) and in order that both sides of 
equation ()5.3|1 are positive we are forced to choose a > 0, 6 < 0. It is remarkable that this is possible 
although J can take either sign, in fact J is not bounded from below, however, it is bounded from 
above by |2a/|6|! 



It follows that 



K 



^pha 



V'd^ L'{J/2) 



L{J/2) 



(5.9) 



e^/a + Jh/2[1 -2K/h]+ P 

e^^[a + Jh/2][l-2K/hf + P 

e^J[a - hV /2\\l - 2K /h] + P 
C 



Vdet(g) 

where the last equality holds when the Hamiltonian constraint is satisfied. 

Using the Brown - Kuchaf key identity V = d/K which holds when the spatial diffeomorphism 
constraint holds and taking the square of ()5.9|) we find 

{^c + (3Y = a + d — (5.10) 



This leads to a quadratic equation for K which we can write as 



{K - A)"^ = ~ B, A 



d 



B 



4a 



(5.11) 



Since the left hand side of ()5.11|) is positive, the right hand side of fl5.11|) is constrained to be positive 
as well. In order to to extend ()5.11|) off the constraint surface we write the two roots of ()5.11|) as^^ 



K = A± X -{A"^ - B + - B\) 



(5.12) 



^Lemma: For real numbers x,y the equation — y implies — {y + \y\)/2. 



Proof: li X — y holds then y > so the second identity follows. The second identity implies the first for y > 0, 
however it implies a; = for y < 0. 
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Since K does not automatically vanish when i? = 0, only the positive sign in (j5.12p is meaningful. 
Moreover, since K is manifestly positive, also the right hand side is constrained to be positive. We 
extend ()5.12j) off the constraint surface by 



A+\j]^{A''-B + \A^-B\ 



A + \j\{A^-B+\A^-B\) 



(5.13) 



which is now manifestly positive. If we want to avoid absolute values then we can write the solution 
explicitly in the form 



K = G(c, d) := - A [(c + pf-a-d] + ^ {- A[(c + /3)2 - a - rf]}^ - ^(6/2)2 (5.14) 

This means in particular that we get the conditions A^ > B, (c+/5)^ > a+d (recall that a > 0, b < 0). 
Hence also (c + /?)^ — {a + d) > 2^Ja\fd, i.e. 

|c + /?|>v^+ya (5.15) 

This condition, if imposed at all x G a, is gauge invariant and invariant under the physical evolution. 
To see this, notice that all quantities in ()5.15p are spatial scalars, therefore ()5.15|) at x is mapped 
to '^{x) where it also holds by assumption. Next, all H{x), x E a Poisson commute with H{M) for 
all M and with H = H{M)m=i in particular. Thus, if H{xY = [det{q)G{c, D)]{x) is positive and 
meaningful for all a; G o" then it will be in every gauge and under the physical time evolution. Thus, 
the extension off the constraint surface ()5.13p which displays H in manifestly positive form is only 
required fo the purpose of quantisation. For classical purposes, ()5.14p is completely sufficient since 
anyway we are only interested in the portion of phase space where fl5.15j) holds. 

Notice that for spatially homogeneous phantom fields the function d is constrained to vanish. 
Hence, for c? = ()5.14j) reduces to 

K=^-^[{c + /3r-a] + \ic + ^r-a\] (5.16) 

which is manifestly non negative and now there are no conditions on c, d, (3, a. We want this expres- 
sion to coincide with for large c in order that the Hamiltonian density equals H{x) = C{x) = 
{^/det{qjc){x) in this limit. This enforces b ~ —2a and \(3\, a small. Specifically we could set f3'^ = a 
in oder to remove the constant term in ()5.16p and b = —2a in order that the coefficient of equals 
unity. We may choose a as small as we want in order to suppress the term linear in c in ()5.16|) . For 
those values the Lagrangean becomes simply L = —y/a[e^/l — I + 5] where 5 = ±1. In case that a 
is not compact it makes sense to choose < a (e.g. (3 = 0) because then ()5.1(j|l vanishes in regions 
where (c + jS)"^ < a. Hence, whatever the fall - off conditions of the fields are at infinity, as long as 
they fall off at all, the support of ()5.16|) will be compact, at least when the diffeomorphism constraint 
holds and when is homogeneous. Therefore, our model can be used, under the assumptions made, 
also if a is not compact and/or has a boundary. For other values of (3, a and in particular if there is 
a cosmological constant present, then the action converges anyway only if a is compact and we then 
also require da = ^ in order to avoid boundary terms. 
Once we have chosen b = —2a, P = 0, ()5.14|) reduces to 



K =^[c^ -a-d]±J{^[c^-a- d]}2 - ad (5.17) 
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It is then tempting to perform the hmit a — so that (j5.17|) becomes 



K =^{[c^ -d]±\c^ -d\) (5.18) 

which has the advantage to be manifestly positive even for d 0. Since for b = —2a = (3 = the 
Lagrangean from which ()5.14|) was derived vanishes identically, ()5.18|) must come from a different 
Lagrangean. Indeed, it comes from the constrained, incompressible dust Lagrangean introduced by 
Brown and Kuchaf in their original work p. Unfortunately, it is not of the form L(//2), rather L = 
— |A(1 — /) where A is a Lagrange multiplier field^^. The canonical formulation of this action leads to 
the contribution to the Hamiltonian constraint given by C^^^^ = |[7r^/ {X^ydet{q)) + {l + V)^ydet{q)X\ 
and elimination of A by its equation of motion results in C^'^^^ = y^7r^(T+l/). One can show explic- 
itly by solving (j5.18|) for c in order to obtain (7?'^'^'^*°™ that (|5.18|) cannot come from a Lagrangean 
of the form L{I/2) by applying the procedure for the inverse problem mentioned at the end of sec- 
tion 13.21 We therefore prefer scalar matter ()5.14|) over the dust matter ()5.18|) which we feel to be 
awkward due to the Lagrange multiplier A. What is awkward about this action is that the Euler - 
Lagrange equations for A require 1 = 1 and do not allow to solve for A while after Legendre transfor- 
mation the constraint equation for A only can be solved by eliminating A^°. and to forget about 1 = 1. 

Let us now consider the other cases. First of all, while a can be as small as we want in or- 
der to suppress unwanted terms in ()5.14|) . the case a = is singular. The case a = 0, that is, 

L(J/2) = ^/Tj2 or L{I/2) = ^JTjf has to be treated independently. For L(//2) = ^JTj2 we 
obtain an expression for the physical Hamiltonian which contains negative powers of c, d which do 

not resemble the usual Hamiltonian. For L(//2) = \JTJ^^ as well as for case iii. in ()5.3p we 
are driven to an algebraic equation of ninth order for K which can no longer be solved by quadra- 
tures, hence this model is ruled out for purely mathematical reasons. Finally, the same analysis 
carried out for case ii. in ()5.3|) leads to K = — c ± \/ (? — d which should be positive, hence c is 
constrained to be negative. Thus C*°* = vr -|- a/ det (g) a/— c ^ \/ (?■ — d and the Hamiltonian would 

be nix] = v/det(g)V-c+ -ci( x) which for small d becomes a/2 det(g)|c|. Thus, the Klein 
Gordon Lagrangen produces the square root of c which is also not what we want. 

Hence the only suitable model corresponds to case i. and we now proceed to explore its proper- 
ties. 

5.2 Physical Properties of the Model 

Let us first check that hi{c,d) := \/G{c, d) of ()5.14|) satisfies the PDE ()3.18|) . Since the calculation 
is not entirely trivial we display here some intermediate steps for the convenience of the reader. To 
simplify the computation we notice that \fG satisfies ()3.18j) with n = 1 ii and only if G satisfies 

^^In their work |Tj they actually used = + WaX"-^ instead of <i> ^ where Wa,X"' are additional six scalar fields 
and defined / = —g^^'^U ^U,y. Then <i> := T is a local clock field and X", a = 1, 2, 3 are local position fields and together 
they are assumed to provide a local coordinate system of M . The functions Wa are also Lagrange multipliers. We 
simply set Wa = here. 

^"If we compare this to the string, then this would be as if we would go from the covariant Polyakov string action to 
the Nabu - Goto string action by eliminating the worldshcet metric (the analog of A here) by its equation of motion. 
However, both string actions are worlsheet covariant while this is not the case here. 
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(pnHl) with n = 2. We find 



'dd 
dG 

dc 



AaR 
h 



2aR 



(G - 6/2) 
G(c + /3) 



where R := VA^ - B, A := b{a + d - {c + /3)2)/4a, B = {h/2fd/a. Then fl^TT^ with n 
h2 := G becomes 



iaR' 



AaR' 



■[{d{G - b/2) - AaRG/b){G - b/2) - G\c + f3f] 

An h 

'■[^G\^[d - (c + (3f] -R) + G{2aR - bd) + {b/2fd] 



'4a 



AaR' 



AaR' 



AaR' 
^AaR' 
^AaR' 
'^AaR' 



An 

'[-jG{A + R){A-R)- aG'^ + G(2ai? - bd) + {bl2fd\ 

An 

'[-^GB - aG^ + G{2aR - bd) + ib/2)^d] 
'l-aG^ + 2aRG + {b/2)^d] 
'[-G^ + 2RG + B] 
'■[GiR+[R-G]) + B] 
'■[iR + A){R- A) + B] 



(5.19) 
(5.20) 

2 and 
(5.21) 



(5.22) 



which is what we wanted to show. 



Next, in order to understand the meaning of e and /3 or 6, we compute the equation of state of 
the model. The energy momentum tensor with our signature convention is given by 



T, 



2 d,/\d^\L{I/2) 



-[g,,{P + ey^^Tlb/2) + '^ ^''^'^ 



2 y/a + bI/2' 



(5.23) 



Energy density and pressure become in the perfect fluid approximation T^^ = prifj^riy +p{g^u + n^ny) 
(with respect to our unit timelike vector field n normal to the foliation introduced in section EI) 



p = T^^n'^n'' = p + e^/{a - bV/2){l - 2K/b) 



P = -ip + g'^T,,) 



l-2K/b 



(5.24) 
(5.25) 



The equation of state "field" is defined by w '■= p/p and becomes for spatially homogeneous for 
which V = 
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where y = a/1 — 2K/h > 1. For /? = we get w = —1 /y'^ i.e —1 < w < independent of e, a. For 
/3 ■ e > we also get —l<w<0. For /5 ■ e < let e = l/^l/i/a- If e > 1 then again — 1 < w < 0. 
For e < 1, w{y) has a maximum at ^/g = 1/e + a/I/c^ — 1 > 1 given by i/i(ye) = (l/z/e)^ < 1 and 
w{l) = —1, w{oo) = hence —l<w<l. Finally for e = 1 we get < w{y) = 1/y < 1. 
Next, the speed of sound, for spatially homogeneous is given by 

a^rt ^ |MZ* = +J_ > (5.27) 

dp op[y)/dy y^ 

independently of /3, ae 

Now if would be the only observable scalar field then we would need p > 0, > for stability 
reasons, hence e = +1 in order that there are no restrictions on the range of K, V and (3 > —\fa 
i.e. e < 1. Then in order to get inflation and the recent accelarated expansion [w < —1/3) of the 
universe (dark energy) we must choose either /? > or at least — i/a < f3. Since, however, the phan- 
tom field is pure gauge and has other purposes, it will not be associated with the physical infiaton 
and/or dark energy. Thus we keep the ranges of e,a,b,P unrestricted up to the requirement that 
a > 0, 6 ~ —2a; (3"^ and a small. We will see however, that the physical evolution equations of the 
next section still impose further restrictions. 

Let us set h = —2a for definiteness. Then we conclude that we obtain the two parameter set of 
Dirac - Born - Infeld type Lagrangeans 

L = y/\det{g)\{-l3 + a^l + g'^-^^^^,,) = a^/\ det{g)\{-l + ^1 + g^-^^^<^^,) + (a - l3)^/\det{g)\ 

(5.28) 

where a := —t^Ja 7^ is a real number with |a| small and (3 is any real number with |/?| small. 
For small I the first term in ()5.28|) becomes to linear order Oa/I d&i{^g^\g^^^ ^^^^ ^^12 which up to the 
constant a is just the usual Lagrangean for a massless Klein Gordon field. Thus, to this order the 
Lagrangean has the correct sign in front of the kinetic term for a < and it becomes a phantom 
in the cosmo logical sense for a > 0. The second term represents a contribution by a — /3 to the 
cosmological constant which is a positive contribution for a — negative. Since the cosmological 
term can always be absorbed into the contribution C to the Hamiltonian constraint by gravity and 
the remaining (physical) matter, a natural choice would be /? = a which would then imply e = 1 
hence < w < 1. 



Summarising, we have found a simple scalar field model which in the physical situation of interest, 
that is, a roughly homogeneous phantom field in order that 0(x) r = const, is a good approxima- 
tion for a physical clock, gives rise to a satisfactory physical Hamiltonian. It decays sufficiently fast 
at spatial infinity (in fact has compact support) when simultaneously the spatial diffeomorphism con- 
straint holds. Decay requirements can be avoided if there is a non - vanishing cosmological constant 
so that a is compact and we impose also c^cr = in this case to avoid boundary terms. 

It would be interesting to improve the model in order that it is applicable also in situations where 
a is bounded and has a boundary. Also, for mathematical reasons in particular in view of quantisation 
(see the next section) it would be more convenient to have a manifestly positive Hamiltonian. A 
possible starting point is the manifestly positive extension (j5.13p off the positivity constraint surface. 
Finally it might be possible to find a model such that H approximates C{N = 1) even when V is 
not close to zero. One way to investigate this would be to solve the inverse problem mentioned at 
the end or section 13.21 We leave this to future research. 
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6 Consequences for Cosmology 



The Hamiltonian constraint for our model is given by 

C"t°t(x) = Ti{x) + ^det{q)G{c, d){x) =: 71 (x) + H{x) and we see that the Hamiltonian becomes 
for small V and large c(a;)^ ^ a^,/?^ approximately H = J^d^xC{x) which is just the integrated 
contribution, of the gravitational and non - phantom like matter degrees of freedom, to the original 
Hamiltonian constraint. It would result from the canonical Hamiltonian by choosing the lapse to equal 
unity, the shift to equal zero and by dropping the phantom field contribution from C*°*, This ex- 

plains why in the presence of our particular phantom field model chosen, evolution with unit lapse and 
zero shift with respect to the canonical, original "Hamiltonian" H^^'^°'^{N, N) = C^°^{N) + D*°*(iV) 
of functions on phase space not involving (p approximately equals the physical evolution of the non - 
phantom like degrees of freedom. 

We will now illustrate the meaning of "approximately" in the context of isotropic and homoge- 
neous minisuperspace models, that is, FRW cosmology. The presence of the phantom field, while 
not directly observable, will still have an important impact on the conceptual (interpretation) and 
technical (matter content) aspects of the FRW equations as well as on their validity. Namely, we 
will see that the FRW equations are only an approximation to the actual physical time evolution of 
observable quantities generated by the physical Hamiltonian. This also serves to explain the formal- 
ism in a simple context without the field theoretic complications^^. 

We consider the experimentally almost confirmed case of a spatially fiat {k = 0) model for which the 
usual FRW line element is given by ds'^ — —dt^ -\- a{tY5abdx"-dx^ where a is called the scale factor. 
Comparing with the general ADM line element ds^ = -[N'^+qabN"'N^]dt^ + 2qabN^dx"-+qahdx"-dx^ we 
read off = 1, A^*^ = 0, qab = aitY^ab ='■ Q^ab- Here we work with dimensionless spatial coordinates 
x"- so that a has dimension cm while the unphysical time (or foliation parameter) t has dimension 
cm. We also take our scalar field to have dimension cm so that / = {d4>/dty is dimensionless. 

We begin by specialising the canonical formulation of GR to isotropic and homogeneous situa- 
tions. The extrinsic curvature Kab — {qab — ^NQ.abj/i'^N) where C denotes the Lie derivative reduces 
to Kab — a da/dtSab, hence the momentum P"^ = ^y^det{q)[q'^q'^'^ — q"-^q'^'^]Kcd conjugate^^ to qab 
becomes P"^ = — d(5"* =: Pq6°'^/3. Here Q,Pq are canonically conjugate. The canonical transfor- 
mation from {Q = a^jPq) to (a, P = ^aPq) reveals P = —6a da/dt. The spatial diffeomorphism 
constraint vanishes identically if spatial homogeneity is assumed and the contribution to the Hamil- 
tonian constraint of the gravitational degrees of freedom and the cosmological term becomes 

C^^^ = i( Vqi^(^[(5-5'«^ _ q-\^)KabK^ - R^'^ + 2A^/dit(^) = + Aa' (6.1) 

where R'^^^ is the curvature scalar of qab which vanishes identically. 

For ordinary matter we will make as usual a perfect fiuid Ansatz for the energy momentum tensor 
Tf_i,u = Pmn^riy + Pm{gfiu + n^riy) with — Sj^ where energy density and pressure are related by 
—Sa'^Pm — d{pmO^) / da. For instance for a Klein - Gordon field Pm — (tt^/o^ + U)/2 where U is its 
potential and thus p^ = (vr^/a^ — U). Since in general c^'^'^^ — —2[G^i, + \Kg^J\n^'ny we find that the 
contribution to the Hamiltonian constraint of gravity and non - phantom matter reads in general 

^^Notice that as usual wc quotient all equations by the infinite coordinate volume / cPx in the /c = 0, — 1 models. 
^^Recall that the gravitational and cosmological action are multiplied by l/(167rG) = 1/2 which explains the factor 
1/2 in front here. 
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C = C^^^^ + a^pm- If there is a phantom present, as we advertised in the present paper then the total 
Hamiltonian constraint reads C*°* = C + aSpPhantom = c + CPhantom_ 

As we showed in section (jS-H) the phantom field contribution to the Hamiltonian constraint is 
given by (remember b = —2a, V = 0) 



^phantom _ ^phantom/^3 _ ^ _ ^^/^^ = ^ _ 6 ^TT^j^.^tom/^' + =: P - eaWT^ (6.2) 

where we now take a > 0, e = ±1 and have introduced the "deviation parameter" 

X := (6.3) 

which will be crucial for what follows. We noticed that the gauge invariant quantity vTphantom is also a 
constant of the physical motion because the physical Hamiltonian H = — vTphantom does not involve 
(it is cyclic) and therefore we have denoted the energy squared constant of motion by ii^^ := TTpj^j^j^^^^jj. 
Since (j6.2p has the same dimension as the cosmo logical constant A we see that a,f3 have dimension 
cm~^, TT, E have dimension cm and x is dimensionless. The phantom pressure is given by 

Pphantom = - TT^ d{a^ ppY,s,ntom) / da = -/5 + ea^tt^/ vW+o^ = -P + £ (6.4) 
3a^ V 1 + X 

This gives an equation of state 

_ Pphantom _ [3 - €0 / y 

"W^phantom — — ^ , l'^-'-'/ 
Pphantom P + £"2/ 



with y = vT+x. The phantom speed of sound is given by 

^2 C^Pphantom/ dy 1 (6 6) 

phantom (ipphantom / C^Z/ 1/^ 

which is always positive independently of (3, a, e. 

The canonical "Hamiltonian" of the theory without the phantom is H'^^^"^ = C . Without the 
phantom, however, it is constrained to vanish and therefore should not be interpreted as a Hamilto- 
nian, rather as a Hamiltonian constraint which generates gauge transformations and not any observ- 
able evolution. With the phantom, C is not constrained to vanish. The physical Hamiltonian, with 
the phantom present then follows from ()5.16|) with h = —2a (remember d = identically in exactly 
homogeneous cosmology) 



H = a'^^[{c + /3)2 -a^] + \{c + (3^ - a^\] (6.7) 

which approaches \C\ = a^\c\ for |c| ^ It vanishes identically when |c-|-/5| < a (this never 

happens on the constraint surface). Hence, as long as c > and c S> |/?|, « the physical Hamiltonian 
H, in presence of the phantom, is in good agreement with the canonical Hamiltonian constraint, in 
absence of the phantom. 

We will now derive the FRW equations from the canonical formalism and compare them with 
the evolution equations of physical observables. We begin with the standard FRW equations which 
consist of a set of two equations. The first one is just the constraint equation C*°* = with P 
eliminated by the equation of motion for a which gives a condition on {da/dt)'^. The second one 
involves d^a/dt^ and is obtained by eliminating dp/dt by its equation of motion. Notice that by 
"equation of motion" we mean actually gauge transformations generated by C*°*. On the other 
hand, we can compute evolution equations generated by the physical Hamiltonian H. 
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1. Gauge Transformations generated by Standard FRW Equations 
The gauge transformation for a is 

^:={C"a} = ^ = 4 (6.8) 

at op 6a 

Inserting ()6.8|) into the constraint equation C*°* = we therefore find the first one of the FRW 
equations 

3( — )Va^ = ^ + Pm + Pphantom (6.9) 

The second equation is obtained by solving ()6.9|) for da/dt, taking the second derivative and 
using the conservation law dp/dt + 3{p + p)da/dt = which gives 

^^Z*^ = ^~^(Pm + ^Pm + Pphantom + 3j9phantom) (6.10) 

This is the second FRW equation for spatial curvature k = and using units for which SttG = 1. 

2. Physical Time Evolution generated by H : Modified FRW equations 

We can now compute the time evolution of physical observables. The general formula ()4.2j] 
specialises to 

^/(^) = E V^H^^ = E Hr = {r- m (6.11) 



n=0 n=0 



with 



H = ^{C + pa^f - a'^a^ (6.12) 



and / can be any function on the cosmological minisuperspace phase space independent of vr, 
because the spatial diffeomorphism constraint vanishes identically. Obviously, just as in the 
full theory (iO/(r) = {if, 0/(r)} = 0{hj}{t). We now see that 

O/(r)^/(r-0) = /(t), t:=r-0 (6.13) 

where 1 1— > /(t), /(O) = / is the solution of the evolution equation for / without the phantom 
when treating the Hamiltonian constraint C as a Hamiltonian. The unphysical time parameter 
t of that unphysical time evolution is now interpreted as t = r — 0. Hence all the cosmological 
evolution equations remain approximately intact (under the restrictions on c made above), 
however, we now have justified why that evolution corresponds to observation and we have 
interpreted the time parameter t as composed of the pure gauge phantom time and the 
physical time parameter r. 

Thus, the phantom has nicely reconciled the mathematical framework (gauge theory) with 
observation (FRW equations). It is pure gauge and one would be tempted to conclude that 
its presence does not have any observational consequences beyond modifying the Hamiltonian 
constraint. This is of course wrong because, at least in our model, we cannot have C = H just 
C ^ H and we now proceed to compute the associated modifications. Hence, what we need to 
do is to repeat the steps (j6.8p - (j6.1(Jj) where the first FRW equation C*°* = C + a^Pphantom = 
must be expressed in terms of observable quantities and instead of the gauge transformation 
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d/dt{.) = .} we now have actual physical evolution d/dT{.) = {H,.}. The physical 

evolution equation for the observable Oa corresponding to a is^^ 



dOa 



{H, Oa} = 0{HM = Oeu.^. = -OshOpHQOo) (6.14) 



Here 



^ = ^— = e ^ =eVlTl/x (6.15) 

Using ()6.14|) in Octot = we find the first modified FRW equation^"^ 

j^(^)' = (^ + ^)l^ + Op™ + Opp.an.oJ (6.16) 

The second again follows with the definition of pressure by taking the second physical time 
derivative of Oa 

^ "^'^^ '1 + 7^)A - + 0,^,_J(1 - A) + 3[0„„ + 0,^,_J(1 + -^))} (6.17) 



^ ^ O" 2 PphantomJV ^ / L PphantomJ V ^ 

Notice that O^ = E'^/{a^Ol). We now interpret ^cTT^ and (IFTTTfl : For O^; ^ oo they look 
exactly like standard FRW equations with an additional phantom matter component. In fact, 
aX Ox = oo these two equations are identical with ()6.9|) and ()6.10|) under the substitution 
(t, a(t)) {T^Oair)). That additional matter component behaves at early times {Oa — > 0) of 
the universe {Ox oo) as 

Pphantom ^ — e-E/a'^, ^phantom — * — (6.18) 

Hence it behaves like dust matter for e = —1 as Oa provided we set /3 = 0. It therefore 
could serve as a dark matter candidate then. At later times, when Oa oo or O^; we get 

Pphantom ^ /3 - 60, Pphantom ^ -/5 + (6.19) 

SO it behaves like a positive/negative cosmological constant for e = — 1 or e = +1 respectively. 
In order to get a positive contribution to the cosmological constant and to retain these inter- 
pretations, we thus should choose /3 = 0, e = — 1 which is exactly of the form used in the model 

However, our purposes are different here: The phantom field is there in order to deparametrise 
the theory and to provide a positive, physical Hamiltonian which approximates C since C is 
used in the standard model with fiat space. Since 0*°* = C + Pphantom = we are forced to 
take Pphantom < in order to have C* > 0. This excludes the e = — 1 case^^ and we must 
use e = +1. Hence, the phantom has negative energy, positive pressure and behaves like 



^■^The convenient observation here is that the map f ^ Of is an automorphism. In particular if / = F{a, P, (p^n) 
for some function F then Of = F{Oa, Op, r, tt). 

^"'Optot = C*°* — foUows because C*°* is already gauge invariant. Here we need to use the generalisation of the 
map f ^ Of mentioned in section 0J 

^^More precisely we could take C"*°* = 7r± -v/C^ — a^a^, i.e. H = ^pVC^ — a^a^ « for small a and a. In order 
that this approximates C we shouUd have C < or C > respectively which requires Pphantom > or Pphantom < 
respectively. However, we should have Pphantom < in order to avoid a big rip singularity as we will show below. Also 
C > in flat space. 
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negative energy dust in the early universe while it apporaches a negative cosmological con- 
stant in the late universe. Since the gravitational contribution to C is negative, we must have 
Pm + Pphantom > in Order that C*°* = 0. Hence the total energy density of matter is positive 
which is sufficient in order to obtain a stable theory, that is, the usual energy conditions are 
satisfied ^2] • In fact the phantom can be compensated for by a k - essence field with positive 
energy |3] or by ordinary positive energy dust matter and a positive cosmological constant term. 

The reasoning here would then be as follows: 

1. Something like a phantom is needed in order to deparametrise the theory and to keep va- 
lidity of standard model physics and the FRW equations etc. 

2. Since the phantom energy is negative we must compensate for it by positive energy matter. 
The simplest way to do this is to add a k - essence field with energy density pk = ^/t^I/cl^ + 7^ 
which just corresponds to an action of the DBI type with a replaced by 7 and e = — 1. The k 
- essence momentum tt^ is also a constant of the motion. 

Thus, both the phantom and k - essence fields are called for by the mathematical formalism. 

Notice that the transition between the regimes where the usual FRW equations retain their 
interpretation as evolution equations of observable quantities is controlled by the deviation 
parameter Ox = E'^ The transition occurs at at O^; = 1 so Oa = \J Eja. By choosing 
a sufficiently small and/or E sufficiently large we can achieve that the transition scale is as 
large as we need. Notice that a is a kinematical parameter of the Lagrangean while E is 
a dynamical constant of motion. However, whatever the value of a, equations ()6.16p differ 
drastically from the standard FRW equations beyond O^; = 1. That is, the universe evolves 
completely differently beyond O^; = 1. The largest corrections at late are of the order of 
1/Oa; oc which are terms normally not considered in the FRW equations. In fact they will 
completely dominate then. 

Let us see what will qualitatively happen at very late times: We will use realistic matter 
composed of dust (baryons), radiation and k - essence (for simplicity without potential term) 
and set 

Pra = B/a^ + R/a^ + ^J-Kl/a^ + ^i"^ (6.20) 

where B,R > are integration constants and 7 > 0. Notice that also vr^ is a constant of 
motion. The physical evolution equations become (we replace Oa by a etc. for the purpose of 
this discussion) 

3{da/dTf/a^ = (1 + 1/X)[A + Pm + Pphantom] (6.21) 

and we clearly need that the right hand side is positive for the entire evolution. During radiation 
domination, i? > is sufficient. During baryon domination we must have B + Ivr^l — \E\ > 0. 
For X — > oo we may want to require A + 7 — a > 0. However, we notice one problem at late 
times: Suppose that the right hand side of ()6.21|) never vanishes. Then da/dr > for all r due 
to continuity. Expanding ()6.2ip around x = we find 

1 22 

Sida/drf/a' = {l + l/x){[A + j-a] + -i^^-a)x + B^V^+Ri^./^Y/' + Oix')} (6.22) 

2 £/^7^ E E 

Even ifA + 7 — a = the right hand side diverges at 2; ^ due to the radiation and baryonic 
terms present. The leading contribution is then given by 

2 

{da/drY = B^a'' =: K^a' (6.23) 
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which can be solved by 

<r) = !- (6.24) 

with K > 0, (5 is a constant of integration and the upper/lower sign corresponds to da/dr > 
0/ < respectively. Hence, if da/dr > for all r then we must take the upper sign and 

6 > and would conclude that the universe reaches infinite size after the finite amount of time 
25/(3k). The evolution would in fact stop there and is called a big rip singularity. This is 
clearly undesirable and the only way to avoid this is to tune the parameters in such a way that 
p can vanish while always being non - negative during the evolution. This is only possible if 
the phantom contribution to the energy budget is negative. The universe would then be able 
to reach maximum size and then would recoUapse. This is granted to happen if the right hand 
side of ()6.21|) kinematically can become negative beyond some critical Xc (dynamically it can 
never happen because [da/drY > 0). In fact it is not difficult to show that we must have 
A + 7 — a < for this to happen. 

One can tune B,R,A,E,iTk,'y,<y such that there is a radiation, dust and positive vacuum 
energy era while still x ^ 1 during which the FRW equations hold. The radiation era holds 
for < a < p , I — ^ =: during which \/x > ^ 1 must hold. The dust era is < a < 

sj B+\-Kk\^ _. during which y/x > ^ 1 must hold. Since A/a < 1 and < a^, both 
conditions are easily satisfied with the observed values for R,B for Itt^I — of the order of B, 
E/B ^ 1 and E/a sufficiently small. Finally, the vacumm era lasts for < a < Oc where ac 
is the value at which ()6.2H) vanishes. In order to see that Xc is smaller than x = 1 at which the 
FRW equations anyway no longer take their standard form we notice that at x = 1 the right 

hand side of ()6.22|) is still larger than A + \J ^"^ + — \/2a which is positive independently 
of the value of A, 7 as long as 7^1/ E"^ > 2. However, if we associate B with baryonic matter 
and Itt^I — E with dark matter then (|7rfc| — E) / B ^ 10. Since B/E <^ 1 we get \Tik\/E ^ 1 
so we should have A + \/ — \/2a > for A, 7, a > subject to A + 7 < a. Fix some 
e > and set y := 7/a and A/a := (1 — e — y). Then we need to find < y < 1 — e such 
that 1 — t — y + a/1 + — -\/2 > 0. Set k := \/2 — 1 + e which is positive and satisfies k <1 
for e < 2 - v^. Then < y = (1 - k^)/(2k) < 1 - e which implies 1 + (1 - e)^ < 2 which 
is identically satisfied for < e < 2. Hence we may choose any < e < 2 — \/2 and then 

7 = a(l — K^)/2fi;, A = a(l — e) —7 where n = a/2 — 1 + e. It is appealing that the cosmological 
constant during the vacuum era A is of the order of a which should be small, thus explaining 
the smallness of the cosmological constant^^ 

Let us summarise once again the observable effect of the phantom: 

1. The physical evolution equations of observable quantities have a standard FRW form for large x 
(small a). 

2. The phantom adds additional matter terms with an equation of state —l<w<0 which evolves 
from — >■ — 1 as a evolves from — *• 00. It therefore acts like in k - essence, just with negative 
energy. However, an additional k - essence field is natural in order to compensate the negative 
phantom energy. 

3. In particular, it is wrong that the physical evolution equations have just the FRW form without 

^^For completeness we also mention a scenario without k - essence matter, that is, tt^ = 7 = 0. In this case the same 
analysis yields the condition < E/ {B — E) ^ 1. Then at a; = 1 the energy density becomes p > A + a{B/E — y/2) 
which is still positive for y/2— 1 — e<A/a< 1 where B — {1 + e)E. Notice that now B — E corresponds to observed 
baryonic matter. 
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phantom matter contribution. It is easy to see that this is due to the fact that 7^ (^ideai ._ _j_ 7^ 
in which case we would have H = C exactly. Rather we have C*°* = C — a/tt^ + a^a^ so that 
C*°* — C^^^^^ = \n\{^yi + 1/x — 1). It is not difficult to see that there is no Lagrangean of the form 
L(//2) which can produce H = C exactly. 

4. As X becomes small, the actual evolution equations differ drastically from the FRW form. The 
transition is roughly at at = ^JeJol where E is the nergy of the universe and a is a parameter of 
the model which can be tuned to be so small that at is way beyond today's value oq. In order that 
the universe has infinite observable life time, the parameters can and must be tuned such that the 
universe in fact recoUapses rather than expanding forever. 

Notice that we are not doubting the validity of Einstein's equations at all. These are completely 
encoded in the fundamental constraint C*°*. There are two of these equations. One is C*°* = 0. 
The other results by computing the gauge transformation da/dt = {C*°*,a}, to solve this for P, to 
compute the gauge transformation dP/dt = {C*°*,P} and to insert this into d'^a/dt^. This results 
exactly in (jfj.lfij) . ()fi.l7j) at x = 00 and with t replaced by r. However, what we critisise is that these 
are interpreted as physical evolution equations. They are not, they are gauge transformations of 
non - gauge invariant, unobservable quantities and not evolution equations with respect to a non - 
vanishing Hamiltonian of observable quantities. What we have done here is to compute the physical 
evolution of observable quantities generated by a physical Hamiltonian. The resulting equations are, 
by a judicious choice of phantom field Lagrangean, in good agreement with the standard equations. 
However we insist that the standard procedure is fundamentally wrong. In particular, the standard 
FRW equations are drastically false in the late universe if our phantom field is realised in nature and 
provides the physical Hamiltonian which generates the evolution that we actually observe. These 
reservations hold of course in full generality in all applications of General Relativity. 

There is another way to look at what is going on here: What we have done in order to obtain 
physical evolution is to use the unphysical gauge transformation d(f)/dt{t) = {C*°*, 0}(^=0(j) and to 
solve the condition (j){t) = r for t. This results in the function t ^ tr which is a non - trivial function 
on phase space. Now we insert the value tr into the unphysical gauge transformation t ^-^ a{t) where 
a{t) solves da/dt{t) = {C^°^ , t} a=a{t) and obtain a{tr). We claim that a{tr) = Oa(r), at least when 
(jtot _ vi^here Oair) solves dOaldrij) = {if, Oa(r)}. To see this, we compute for any function / 
with equation of motion df /dt = {C*°*, /} 

{c'^\f{u)} = {c'°\f}f=m,t=u + {df/dt)t=tAc'°\u} = {df/dt)t=tAi + (6.25) 

Now choose / = in fl6.25|) and use that r = is a constant function. Then use / = a in fl6.25|) 
to see that a(tT-) is an observable. Now on the constraint surface 

= {C"*°\/(t.)} = {7iJ}f^t),t=u + {HJ}m,t=t^ + {df/dt)t=tA{ix.U} + {H,U}] (6.26) 

Choose / = and use that H does not depend on vr, then 

= 1 + {d<P/dt)t=t. [{tt, t,} + {H, U}] (6.27) 

Insert ()6.27p into ()6.26|) with / = a to obtain 

= {H,a}a=a{t),t=tr - {da/dt)t=tj[{d(t)/dt)t=t,] = {H,a}a=a(t^) - da{tr)/dT (6.28) 

Thus a(tr) and Oa(r) differ at most by a constant. 
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What happens here is that t CL{t), t ^ (j){t) is the parametrisation of a trajectory in phase 
space (here restricted to the a, plane) which is obtained by solving the equation (pit) = r for t 
and inserting this into a{t) so that we arrive at r i— *• a{tr)- We have deparametrised the description 
and now are dealing with the only physically meaningful object, the trajectory itself and not some 
random parametrisation thereof. Any reparametrisation t = t{t') with dt/dt' > 0, that is, a gauge 
transformation, changes the unphysical functions r{t), (j){t) but results in the same physical trajectory. 

What consequences does this have for the FRW line element 

ds^ = -dt^ + a{tYdx^dx''5ab (6.29) 

which is also expressed in terms of the unphysical quantities t, a{t)l Let us express the line element 
in terms of r by applying the diffeomorphism t := tr- In these coordinates ()6.29|) becomes 

ds'' = -rfr2[^(t)]2^,^ + a{uYdx''dx''5ab = -dr^l + ^^^) + aiuYdx'^dx'Sab (6.30) 

which is again a FRW line element, now expressed in terms of observable quantities, at least for small 
a. For large a, ()6.30|1 is no longer of FRW form. Again the deviation parameter x has appeared and 
shows that for sufficiently small a the metric ()(j.3()|l expressed in terms of observable quantities is 
well approximated by the usual FRW form even today. 



7 Conclusions and Outlook 

It has been known for a long time that the problem of time can be solved in principle by the relational 
framework due to Rovelli and others. This has never been appreciated as much as it should have 
been because, while the conceptual, physical framework was clear, the analytical implementation 
remained largely undeveloped for a long time. With the appearance of [S], analytical methods became 
available for the first time. Still the framework, in its full generality, remains discouragingly difficult 
in particular when applied to General Relativity due to the complexity of the analytical expressions 
which involve summing an infinite number of infinite series, an inversion of infinite dimensional 
matrices and the computation of an infinite number of different, iterated Poisson brackets^^. 

The first main message of the present paper is that when adding appropiate, albeit hypothetical, 
matter, the complexity of these formulas is drastically reduced. In contrast to the general case, 
there is only one series to sum, there are no matrices to invert, there is only one kind of iterated 
Poisson bracket to compute. Hence the formulas that we obtain are remarkably simple. In fact, the 
classical time evolution in a background dependent theory, say in QCD on Minkowski space, of some 
observable O such as a Wilson loop function, would also be given by the series 

oo „ 

^(^) = E^{^QCd,0}h (7.1) 

n=0 

where -ffqco is the QCD Hamiltonian. Comparing with (j4.2|) we see that the complexity of the 
classical time evolution in both theories is comparable! 

^^Drastic simplifications occur as far as the number of relevant constraints is concerned (four instead of infinitely 
many) when reformulating GR in terms of coordinates that are spacetime scalars ^3] . While the resulting observables 
are relatively simple (although still inversions of non - trivial matrices take place) and physically intuitive, the field 
variables that one uses are complicated compounds of the the canonical fields and the observables involve polynomials 
of those evaluated at one specific spatial point. In quantum theory, these observables are therefore presumably too 
singular because they involve the product of operator valued distributions evaluated at the same point. 
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The second main message is that, in contrast to the general case, the physical observables we 
obtain are strong observables and there is just one natural, physical Hamiltonian which does not 
depend on the physical time parameter. That Hamiltonian is (constrained to be) positive and at 
least in the physically interesting region in phase space, that Hamiltonian reduces to the canonical 
Hamiltonian that one usually uses in cosmology and the standard model when the metric is fiat. In 
fact, we manage to completely deparametrise the system irrespective of the other matter present. 

The third main message is that the scalar type of matter that we considered here, from the 
mathematical (to be able to solve algebraic equations) and physical (to obtain a physical Hamiltonian 
which is close to that of the standard model) perspective naturally leads to Dirac - Born - Infeld 
(DBI) negative energy phantom fields with constant potential. This negative energy phantom must 
be compensated for by positive energy matter, most naturally by a k; - essence field. Such matter 
was discussed independently in the cosmology literature in order to provide a candidate for inflation 
and dark energy. Hence the scalar matter we consider here might actually really exist! 

The fourth main message is that, at least for the scalar model we have used here and for which 
we gave strong motivations, the usual interpretation of the cosmological framework, although funda- 
mentally wrong because gauge transformations of gauge dependent objects are interpreted as actual 
physical evolution equations of observables, remains valid when analysed in the correct way, that 
is, by computing the physical evolution of gauge invariant observables. The domain of validity of 
these equations can be tuned to be arbitrarily long, however, it is manifestly finite when using the 
physical time parameter corresponding to the physical Hamiltonian. The actual evolution at late 
times apparently leads to a recollapse. 

The future application of this framework lies of course in the quantum theory. The framework 
presented here, as well as any other application of the relational Ansatz so far, is purely classical. 
In order to promote the framework to the quantum theory, the functions H should be promoted to 
positive self - adoint operators and the functions and H{M) to self -adjoint operators^^. If we 
find operator orderings such that [H{M),H{M')] = 0, [H{M), Ht] = ihH{M) then the quantum 
observables are given by 



where the unitary operators displayed are defined by the spectral theorem. They manifestly com- 
mute, under the assumptions made, with exp{iH[M)/h) and the operator ordering problem for the 
observables would be solved^^ 

Given these assumptions, the fact that then a positive, fundamental Hamiltonian is available 
could enable one to solve the vacuum problem in quantum cosmology: Namely, in usual semiclassical 
quantum cosmology one neglects quantum gravity and applies the framework of quantum field theory 
on curved spacetimes ^HI- The issue is that in cosmology the background metric is not stationary 
and therefore the problem becomes, roughly speaking, to choose a point of unphysical time and at 
that time a definition of annihilation operators (for the free fields) suggested by the energy density 
function of the matter in question in order to select a vacuum state. This is highly ambiguous 
and the cosmological evolution does not keep the vacuum intact but rather causes constant particle 

^^These operators are supposed to be spatially difFeomorphism invariant, see |10| for a quantum implementation of 
the difFeomorphism group within LQG. 

^^If one cannot find a model in which all expressions of which one has to take the square root are manifestly positive, 
then we may be able to compute the spectrum of the (regulated) operators without the square root and restrict the 
Hilbert space to the "subspace" on which all of them take positive (generalised) eigenvalues. This is possible because 
the operators are supposed to commute. The square root would then be well defined on that subspace as has been 
pointed out in pp. Alternatively one can try to use the manifestly positive substitute expressions H5.13|l . 
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production. On the other hand, if one has a fundamental Hamiltonian at one's disposal, then it is 
natural to define a vacuum state as a minimal (zero) energy (eigen)state. This would circumvent 
this problem of initial conditions. 

Having a physical Hamiltonian and physical observables at one's disposal one can also hope to 
develop physical scattering and S - Matrix theory. Namely, while we drastically simplified the rela- 
tional framework, it will be still very hard to compute 0/(r) explicitly to all orders. Here one will 
use the series in order to perform perturbation theory in the way outlined in jTT], say within the 
framework of LQG: Given approximate physical states which can be obtained by using semiclassical 
techniques of LQG ^Hj, we can concentrate them on regions in phase space where (f){x) ~ r. Then 
expectation value computations of physical observables can be terminated after a few terms in the 
power series and only a small number of iterated commutators has to be computed. This should 
work especially nice in applications to quantum cosmology [II] within LQG. 

In summary, there is much left to do in order to make this framework practically applicable and 
it is worthwhile to explore the space of Lagrangeans which lead to deparametrisation further. How- 
ever, we feel that conceptually the framework is quite clear, the complexity has been drastically 
reduced, its validity has been checked in a cosmological setting and the remaining technical tasks to 
be solved have been identified. 
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